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ABSTRACT 

In this paper, we present spectre-photometric models for galaxies of different mor- 
phological type whose spectral energy distributions (SEDs) take into account the 
effect of dust in absorbing UV-optical light and re-emitting it in the infrared (IR) . 
The models contain three main components: (i) the diffuse interstellar medium 
(ISM) composed of gas and dust whose emissio n and extinction properties have 
already been studied in detail bv iPiovan et al.l l)2006|) . (ii) the large complexes 
of molecular clouds (MCs) in which new stars are formed and (iii) the stars of 
any age and chemical composition. 

The galaxy models stand on a robust model of chemical evolution that as- 
suming a suitable prescription for gas infall, initial mass function, star formation 
rate and stellar ejecta provides the total amounts of gas and stars present at any 
age together with their chemical history. The chemical models are taylored in 
such a way to match the gross properties of galaxies of different morphological 
type. In order to describe the interaction between stars and ISM in building up 
the total SED of a galaxy, one has to know the spatial distribution of gas and 
stars. This is made adopting a simple geometrical model for each type of galaxy. 
The total gas and star mass provided by the chemical model are distributed over 
the whole volume by means of suitable density profiles, one for each component 
and depending on the galaxy type (spheroidal, disk and disk plus bulge). The 
galaxy is then split in suitable volume elements to each of which the appropriate 
amounts of stars, MCs and ISM are assigned. Each elemental volume bin is at 
the same time source of radiation from the stars inside and absorber and emitter 
of radiation from and to all other volume bins and the ISM in between. They 
are the elemental seeds to calculate the total SED. 

Using the results for th e properties of the IS M and the Single Stellar Popu- 
lations (SSPs) presented bv IPiovan et alJ l)2006D we derive the SEDs of galaxies 
of different morphological type. First the technical details of the method are 
described and the basic relations driving the interaction between the physical 
components of the galaxy are presented. Second, the main parameters are exam- 
ined and their effects on the SED of three prototype galaxies (a disk, an elliptical 
and a starburster) are highlighted. The final part of the paper is devoted to assess 
the ability of our galaxy models in reproducing the SEDs of a few real galaxies 
of the Local Universe. 

Key words: Galaxies - dust, extinction, emission; ISM - radiative transfer; 
Galaxies - formation, evolution, population synthesis; Galaxies - ellipticals, spi- 
rals, starbursters 



1 INTRODUCTION 

In the early thirties of the past century, nearly at the 
same time of the Hubble discoveries (external galaxies 
and expansion of the Universe), another important dis- 



covery albeit less spectacular was made by Trumpler, i.e. 
the existence of the interstellar dust, whose implications 
could not be even imagined seventy years ago. He sug- 
gested the existence of a solid component, made by par- 
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tides of many irregular shapes and dimensions, mixed 
and associated to the ISM and absorbing the stellar ra- 
diation: the interstellar dust. It was soon clear that the 
dust grains diffusing and absorbing the light emitted by 
the stars, in particular in the UV-optical region, cannot 
be neglected when measuring the light emitted by distant 
stars and galaxies. 

In the local universe a turning point was the IRAS 
survey that discovered ten of thousands of galaxies, the 
major part of them too weak to be included in the optical 
catalogs, emitting more energy in the infrared than in all 
the other spectral regions. Galaxies with IR luminosity 
as high as Lib. > lO 11 !/©, are the main population of 
extragalactic objects in the Local Universe. 

Which is the reason of this huge IR luminosity? 
What powers the IR emission? Observations show that a 
great deal of the IR luminosity is emitted by dust, which 
absorbs the UV-optical light emitted by stars formed in- 
side the MCs in huge bursts of activity and re-emits it in 
the IR. About the 30% or more of the stellar light in the 
Local Universe is reproc essed by dust in the IR-submm 
range l|Soifer et al.lll987h . 

An important new element in the study of high- 
redshift universe has been the detection by the COBE 
satellite of a FIR and sub-mm background radia- 
tion of extragalactic origin, the so-called Cosmic In- 
fraRed Background (CIRB). This radiation implies 
that galaxies in the past should have been much 
more active in the far-IR than in the optical, and 
very luminous. It is likely that dust plays a fun- 
damental ro le in shaping the SEP of these galaxies 
1996; Guiderdrari^^AJ^ffi]; iFixsen et alJ 
Il998l: iPozzetti et al.lll99gf). Observations 



Sg^^il 
Hausci 



with SCU BA iBareer et al]|l999l: iHughes et aljll99Sfr . 
ISOP HOT llPueet et al. ll99fll:lDole et alJll999h and ISO- 



CAM fclbaz et alJll99SUl999l) have measured the CIRB 
at selected wavelengths trying to detect and identify the 
sources of this radiation. Even if it is not easy to identify 
these objects and to measure their redshift jLillv et al] 
Il99fll : iBarger et, aljll99flt) . from the observations it turns 
out that the IR luminosity seems to be emitted be the 
counterparts at high redshift of the local LIRGs and 
ULIRGs, (ongoing high star formation, obscuration and 
reprocessing of the stellar radiation by a dusty environ- 
ment). 

Observations o f high-z galaxies at z = 3 — 4 
JSteidel et al] Il996t iMadau et all Il996l) confirm that 
these objects are characterized by strong obscuration and 
emission by dust in the IR, in such a way that only taking 
into account their emission all over the spectral range it 
is possible to study their properties. Going to even higher 
redshifts, dust still plays a fundamental role even in ob- 
jects at redshift z > 4, as indicated by the observations of 
QSO s and galaxies (iMcMahon et al]ll994t lOmont et all 
ll99rllSoifer et al.lll99Sh . 

It clearly appears how it is no longer possible to leave 
dust aside in studies of the Milky Way, the Local and the 
High Redshift Universe. 

Stars and dust are therefore tightly interwoven, even 
if the presence of dust is more spectacular in galaxies 
with strong star formation. In disk galaxies, with active, 
mild star formation, dust is partly associated to the dif- 



fuse ISM, partly to the molecular regions with star for- 
mation, and, finally, partly to the circumstellar envelopes 
of AGB stars. The contribution from all the three kinds 
of source is evenly balanced. In starburst galaxies, the 
situation is the same as above but now the key role is 
played by the regions of intense star formation. In el- 
liptical galaxies, which show weak emission in the FIR 
(IRAS), dust is essentially associated to AGB stars of 
small mass that continuously loose matter refueling the 
ISM of gas and dust. The great importance of dust in 
relation to the formation and evolution of galaxies is evi- 
dent. Dust strongly affects the observed SED of a galaxy, 
hampering its interpretation in term of the fundamental 
physical parameters, such as age, metallicity, initial mass 
function (IMF), mix of stellar populations, star formation 
history (SFH). Determinations of luminosity functions, 
number counts and many large-scale relations, such as 
the Tully-Fisher just to ment ion one, are als o affected by 
the dust content of galaxies dOala£ttjl2QQl|) - The st rong 
effect of dust, both for l ocal teell fc Kennicutll200lh and 
high redshift galaxies dMadau ^t^lT^9ttlSte*idel et al] 
1999), is the reason why the evolution with the redshift 
of the galaxy sta r formation rate (SFR) is still a matter 
of vivid debate jMadau et al.1 11999: ISteidel et aT1ll999l: 
IBarger et aljl200rft. 

To get precious informations on the process of galax- 
ies formation and evolution, we need to measure the SFR 
in galaxies at different redshifts, to understand when and 
how galaxies form their stars. To this purpose, many 
wide-field and all-sky surveys are currently running or 
have been just co mpleted, e.g. the G alaxy Evolution Ex- 
plorer [GALEX] llMartin et al.lll99^ and the Sloan Dig- 
ital Sky Survey fSDSSl jYork et al] l2000fr in the UV- 
optical range; the Tw o Micron All-Sky Survey [2MASS] 
ijSkrutskie et aljll997l) and the De ep Near-Infrared Sur- 
vey of the Southern Sky [DENIS] fepchtein et aljll997l) 
in the near infrared. Recently, the Spitzer Space Tele- 
scope has open the gate to dedicated studies aimed at im- 
proving our knowledge of the Universe in the middle and 
far infrared (MIR and FIR), e.g. the SI RTF Wide-Area 
Infra red Extragalactic Survey (SWIRE - lLonsdale et alJ 
120031) and the SIRTF Nearby Galaxies Survey (SINGS 
- iKennicutt et, alll2003l) . The IR data will increase even 
more with the coming ASTRO-F mission dPearson et al] 
2004) and the advent of the James Webb Space Telescope 
(JWST). Combined with other astronomical databases, 
they provide a huge amount of UV-optical and IR data for 
millions of galaxies. Both spectral intervals are strongly 
affected by dust. 

To fully exploit the information on the physical prop- 
erties of the observed objects, this wealth of data must be 
accompanied by the continuous upgrade of the theoreti- 
cal models of galaxy formation and evolution. This means 
that the effects of dust must be included in the SEDs 
of single stellar populations (SSP), assemblies of coeval 
stars with the same initial chemical composition, and in 
the SEDs of galaxies, composite systems made of stars of 
any age and chemical composition, gas and dust. Current 
theoretical models for both SSPs and galaxies in presence 
of dust and still unsatisfactory for many different reasons: 
the lack of theoretical spectra for very cool likely dust- 
rich stars, the poor knowledge of the dust properties in 



© 2005 RAS, MNRAS 000, 000-000 



Galaxies spectra with interstellar dust 3 



high metallicity environments and of the production of 
dust itself, the relation between gas and dust content. 

With a few exceptio n s (e.g iBres san et all |l998l: 
iMouhcine fc Lanconl 120021 : IPiovan et alJ 1200.1 l200ft l. 
the light emitted by SSPs is modelled neglecting the 
presence of dust around their stars , both for AGB and 
voung stars fe.g. |Be rtelli et alJll99l iTantalo i et al.lll99Sl : 
lOirardi et all 120021: iBruzual fc Charlotll2003ll When we 
fold many SSPs using the classical evolutionary popula- 
tion synthesis technique (EPS) we simply convolve their 
fluxes with the SFH of the galaxy. Many spectropho- 
tometric models of galaxies are built in this way: there 
is no dust at the level of SSPs and again no dust at 



the level of the galaxy model (see e.g. |Arimoto & Yoshii 


1987 


: lArimoto & Tarrabl ll99fJ: IBruzual & Chariot 


1995 
1999 


; iLeitherer & Heckman Il995t iLeitherer et al 


; ITantalo et alJ 1999: Kodama & Arimoto 


1997]; 


Kodama et al.lll998l: iTantalo et al.lll998!: iBuzzoni 


2002. 


200511. There are many reasons for this lack of realistic 



dusty models. First, before IRAS (1983) and COBE 
(1989), the important role played by dust in the galaxy 
SEDs was not fully appreciated and consequently only 
the stellar component was take n into conside ration 
llTinslevI Il972t ISearle et alJ Il973l : iHuchral Il977fl . Sec- 
ond, the inclusion of the dusty component and the 
IR emission, in particular, require a higher level of 
sophistication of the models. Indeed one has to develop 
a 3D-model in which the sources of radiation and the 
emitting/absorbing medium are distributed; one has to 
face and solve the problem of the radiative transfer; 
one has to know optical properties of the dust; one has 
to simulate in a realistic way the interactions among 
the various physical components of a galaxy and the 
computational cost is often very high. 

Despite the above difficulties, many efforts have 
been made to develop more and more complex and 
refined models, trying to take into account both the 
effects of attenuation and emission by dust. Considering 
the complexity of the problem, many even sophisticated 
studies are limited to the UV-optical region of the 
spectrum and consider or suggest and discuss only the 
attenuation of the ste llar radiation by dust at various 
levels of detail (IGuiderdon i fc Roc ca-Volmerangel 
1981: IBressan et al.l Il994t ITantalo et alJ I l996|: 
Fioc fc Rocca-Volmerangel Il997l: ITantalo et all Il998l ; 
Buzzonill200l IBruzual fc Charlodfeooa iBuzzonihoOEft . 
In some cases, the emission of dust in the IR/sub-mm 
range is considered iRowan-Robinson fc Crawfordlll989l : 
IGuiderdoni et all Il998t e.g.l but no detailed model 
of the stellar sources whose radiation is reprocessed 
by dust is developed. In this case there is not a clear 
relationship between the sources of UV flux and the 
reprocessed IR flux. There are also models that include 
at different level of sophistication the effects of obscura- 
tion and emission by dust, but many of them have been 
developed to study particular objects or aspects of the 
radiative transfer, or the effects of the space distribu- 
tion of the dust fefsta thiou fc Rowan-Robins onl fl 995 ; 
Bianchi et all Il996t IWise fc Silval Il996t ICimatti et all 
19971: iGordon et alJ Il997t iFerrara et all Il999l: 
Gordon et alJ l200lt iMisselt et alJ l200ll iDopita et alJ 
20051 : iDopitaf [20051 They cannot be applied to a gen- 



eral spectrophotometric study of galaxies of different 
mor phological type. The most recent models are those 
bv [ D evriendt et al] Jl99Sh. IDevriendt fc Guiderdonil 
|2000[) ISilyaet alJ fl998ft . ISilval Jl999l). ITakagi et al 
12003 . 12004 120031) The model of ISilva et alJ <ll99Sft has 
been later u pdated and improve d by including the ra- 
dio emission (IBressan et al.N2002ll the nebular emission 
jPanuzzo et all2003Tr and the X-ray emission iSilva et alJ 
|2003j). They present some important differences in the 
way they approach t he pr oblem. For instan c e, the 
models of ISilva et alJ (Il998l) and ITakagi et all (120031) 
are merely theoretical and all properties are derived 
from a few imp ortant assumptions and/or ingr edients, 
whereas those of IDevriendt fc Guiderdonil i200(J) rescale 
the SEDs to match the average IRAS colors. 

As in the meantime, much progress has been made 
in many aspects of this complex problem, for instance 
bette r understanding of t he dust properties (iLi fc Draind 
l200ll : iDraine fc y|200lh . more complet e grids of stel- 
lar evolutionary models a nd isochrones iSalasnich et alJ 
120001 : iGirardi et aljEoOol). new libraries of stellar spec- 
tra at high r esolutio n ijZwittCTet al Jl2004T) . better chemi- 
cal models iPort inari et alJll998f) . new SSPs with dust 
jPiovan et alJ 120031 l200rj) . an d fina l ly mo re detailed, 
multi-phase chemical models jDwekl Il998f) , we intend 
to present here a new model of population synthesis in 
dusty conditions trying to take advantage of some of 
th ese advancements . It follows the theoretical approach 
of ISilva et all 1^998), but yet improves upon important 
aspects. Major changes are related to the chemistry and 
optical properties of the dust (going from using updated 
cross sections for absorption and emission processes, to 
including new grain size distribution laws, to more accu- 
rate description of the radiative transfer for star forming 
regions), to better models of star formation and chemical 
enrichment histories in galaxies including the most recent 
chemical yields and the effect of galactic winds (whenever 
required), and finally to more recent libraries of SSPs to 
calculate the photometric properties of the galaxies. The 
resulting galaxy SEDs go from the far UV to the far IR. 



2 STRATEGY OF THE STUDY 

The model we have adopted is shortly summarized in 
Sect. |3] where first we define the galaxy components we 
are dealing with, i.e. bare stars, stars embedded in MC 
complexes, and diffuse ISM ("Sect. |3~TI ; second we out- 
line the recipes and basic equations for the gas infall, 
chemical evolution, initial mass function and SFR (Sect. 
l3~2"ll ; third we describe how the total amounts of stars, 
MCs and ISM present in the galaxy at a certain age 
are distributed over the galaxy volume by means of suit- 
able density profiles, one for each component (Sect. l3~3l 
that depend on the galaxy type: disk galaxies (Sect. 
13. 3. H . spheroidal galaxies (Sect. 13. 3. 211 . and compos- 
ite galaxies with both disk and bulge (Sect. 13.3.31 . In 
Sect. 0] we briefly summarize some use ful basic relation- 
ships/equations bv lPiovan et all J2OO6) on modelling the 
dusty ISM. In Sect.[^]we explain how the SEDs of galax- 
ies of different morphological type are calculated. First 
the technical details of the method are described and 
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basic relationships/equations describing the interaction 
between the physical components of the galaxy are pre- 
sent ed. Second, we shortly de scribed the SSPs library in 
use l|Piovan et alJ I2003L l200ff) and the space of the pa- 
rameters. Third, three prototype galaxies (a disk galaxy, 
an elliptical galaxy and a starburst galaxy) are used to 
show the effects on the galaxy SED of varying the pa- 
rameters. This is presented in Sect. 15. 2H for an elliptical 
galaxy, in Sect. l5,2~2"l for a disk galaxy, and finally IS". 2, 31 
for a starburst galaxy. The final part of the paper is de- 
voted to assess the ability of our model in reproducing the 
SED of a few galaxies of different morphological type be- 
longing to the Local Universe. In Sect.HJ we present two 
late-type spiral galaxies, in Sect. |7| we show two early- 
type galaxies and finally in Sect. [S] we examine two well 
studied starburst galaxies. Some concluding remarks are 
drawn in Sect. |5] Finally, the mathematical derivation 
of some expressions presented in this paper are given in 
Appendices A, B, and C. 



3 GALAXY MODELS 

3.1 The main components and outline of the 
model 

First, the galaxy models we are dealing with contain at 
least three components: 

1) The diffuse ISM, made of gas and dust. The physical 
properties of the ISM with dust have been thoroughly 
discussed bv IPiovan et alJ i2006fl . In Sect.|lJ we briefly 
summarize the main properties and useful equations of 
the model we eventually choose for the ISM of galaxies. 

2) The large complexes of MCs in which active star 
formation occurs. In our model we do not take HII regions 
and nebular emission into account. Very young stars are 
totally or partially embedded in these parental MCs and 
the SEDs emitted by these stars are severely affected by 
the dusty MCs en vironment around them and skewed 
toward the IR. fsee lPiovan et al.ll200a . for more details) 

3) The populations of stars that are no longer embed- 
ded in their parental MCs. These stars can be subdivided 
in the intermediate age AGB stars (from about O.f Gyr 
to a few Gyr) that are intri nsically obscu r ed by their own 
dust shells as described in lPiovan e t al. (2003f), and the 
old stars which shine as bare objects. 

Second, the amount of stars and gas (and its compo- 
nents) present in a galaxy at any age must be the result of 
the star formation activity framed in a suitable scenario. 

Third, we have to adopt a suitable scheme for the 
3D distribution of the three components in the galaxy 
volume, in order to describe their interaction, and to cal- 
culate the transfer of the radiation across the ISM of the 
galaxy. 



3.2 Gross history of star formation and 
chemical enrichment 

The star formation and chemical enrichment histories of 
the model galaxies are desc r ibed b y the so-called inf all- 
model developed by IChiosil (Il980t) , and ever since used 



by many authors among wh o m we recall iBressan et ail 
Jl994).lTantalo et alJ (ll99eT ). iTantalo et alJ (ll99ST> . and 
|PcTtmari^t^Ln^9a) . In brief, within a halo of Dark 
Matter of mass Mo, radius Rd, and hence known gravi- 
tational potential, the mass of the luminous matter, Ml , 
is supposed to grow with time by infall of primordial gas 
according to the law 



dM L (t) 
dt 



Mq exp 



(1) 



where r is the infall time scale. The constant Mo is fixed 
by assuming that at the present age to the mass Ml (t) 
is equal to Ml (tG), the luminous asymptotic mass of 
the galaxy (see also for more details ITantalo et alJ il996t 
1998?)): 



Mo 



Ml (to) 



t [1 - e(-*G/T)] 
Therefore, the time variation of the baryonic mass is 



(2) 



M L (t) 



M L (t c 



r[l 



;/r)] 



(3) 



Indicating with M g (t) the mass of gas at the time t, 
the corresponding gas fraction is G (t) = 



M 9 (t) 



Denot- 



M L (t G ) ■ 

ing with Xi (t) the mass abundance of the i-th chemical 
species, we may write Gi (t) = Xi (t) G (t) where by defi- 
nition Xi — 1. The general set of equations governing 
the time variation of the generic elemental species i in 
presence of gas infall, star formation, and gas restitu- 
tion by d ying stars has been intro duced and numerically 
solved by iTalbotfc Arnettl dl97lT) for closed models and 
modified bv IChiosil (|l980j) to include the infall term for 
open models: 



dGj (t) 
dt 



= —Xi (i) * (i) + 



(4) 



+ 



J Ml 



dt 



Gt (t) 



where * (f) is the SFR, $ (M) is the IMF (by mass), Mi 
and M u are respectively the lower and upper bounds of 
the IMF, tm and t' M = t — tm are the lifetime and the 
birth time of a star of mass M. Ri(t' M ) is the fraction of 
a star of initial mass M that is ejected back in form of 
the species i. The first term on right-end side represents 
the depletion of the species i from the ISM due to the 
start formation; the second term represents the growth 
of the species i ejected back to the ISM by stars. The last 
term [-^Gi (t)] is the gas accretion rate by infall. In 
the above formulation only the ejecta of single stars are 
included. 

It has been modified to include the contribution 
of Type la supernovae assumin g that the or i ginate 



in close binary sys tems jMatteucci fc Greggiol 119861 
iPortinari et al.Nl99Sft which are supposed to obey the 



same IMF of single stars: 
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dd (t) 



dt 

Mr 



-Xi (t) * (t) + 



(5) 



y(t' M )Ri(t' M )<S>(M)dM- 



+ (I -A) 9(t' M )Ri(t' M )0(M)dM + 

J M nl 



f(t')Ri(t M )<S>(M)dM- 



Mbv 



A I -4^/ F(n)*(t Ml )Ri(t' Ml )dndM B + 



EsNH ■ A 



M B jo 

mbu $(m s ) r - 5 

i Bl M B J m 



dt 



Gi (t) 



Mbi and Mg„ are the lower and upper limit for the mass 
of the binary system (Mb is the total mass of the bi- 
nary system, Mi and M2 represents the mass of the pri- 
mary and secondary star), F (fj,) is the distribution of 
the fractionary mass of the secondary p — M2/MB, pi 
is the minimum value of this mass ratio and Esnh are 
the ej ecta of SNae la, assum ed to be independent of Mb 
or p iPortinari et al.l|l998l) . In the range between Mbi 
and M bu of binary systems suitable to become a SN la, 
the contribution of single stars (the fraction 1 — A of the 
total) is separated from that of binaries producing SNas 
la (the fraction A of the total). 

The first three integrals on the right-end side rep- 
resent the contribution of the ejecta of single stars. The 
fourth integral represents the contribution of the primary 
star in a binary system, assumed to be unaffected by its 
companion, as far as the released chemical ejecta are con- 
cerned. The fifth term is the contribution of SNas la ex- 
ploding when the secondary star pours all its ejecta on the 
primary star. The factor multiplying Esnh is Rsni (t), 
the ra te of SNae la as described in [ Greggio fc Renzinil 
(I1983D . Following lPortinari et aD (ll99S^ we have adopted 
Mbi = 3M Q , M B u = 12M© and A = 0.2. The distribu- 
tion function of the fractionary mass o f the secondary is 
F (p) = 24p 2 (iGreggio fc Renzinilll983T) . 

The SFR, i.e. the number of stars of mass M born 
in the time interval dt and mass interval dM, is given 
by dN/dt = \V (t) $ (M) dM The rate of star formation 
\]> (t) is the ISchmidtl lll959l) law adapted to our model 
W (t) = vM g (t) h which, normalized to Ml (^g), becomes 



(6) 



The parameters v and k are extremely important: k yields 
the dependence of the star formation rate on the gas con- 
tent. Current values are k = 1 or k — 2. The factor v 
measures the efficiency of the star formation process. 
In this type of model, because of the competition between 
gas infall, gas consumption by star formation, and gas 
ejection by dying stars, the SFR starts very low, grows 
to a maximum and then declines. The time scale r (eqn. 
roughly corresponds to the age at which the star for- 
mation activity reaches the peak value. 

The chemical models are meant to provide the mass 



of stars, M* (i), the mass of gas M g (t) and the metallicity 
Z (t) to be used as entries for the population synthesis 
code. 

We also introduce in the model composite galaxies 
made of a disk and a bulge. In this case the mass of the 
galaxy is the sum of the two components. The disk and 
the bulge are assumed to evolve independently and for 
each component the evolution of its M f (t), M g (t) and 
Z it) will be followed. 



3.3 The spatial distribution of stars and ISM 

In the classical EPS, the structure of a galaxy has no 
effect on the total SED which is simply obtained by con- 
volving the SSP spectra with the SFH. The galaxy struc- 
ture is mimicked by considering different SFHs for the 
various morphological types and/or comp onenents of a 
galaxy, e.g. bulge and disk (|Buzzon3[2002). This simple 
convolution is no longer suitable to be used when the ISM 
and the absorption and IR/sub-mm emission by dust are 
taken into account. In this case the spatial distribution of 
the ISM, dust and stars in the galaxy should be specified. 
To this aim we start from the observational evidence that 
the spatial distribution of stars and ISM depends on the 
galaxy type. In our models we consider three prototype 
distributions, i.e. a pure disk, a spheroid and combination 
of the two to simulate late- type (with no bulge), early- 
type (classical ellipticals), and intermediate-type (with a 
prominent bulge) galaxies, respectively. 



3.3.1 Disk galaxies 

The mass density distribution of stars, p*, diffuse ISM, 
pM, and MCs, pMC, inside a disk galaxy can be repre- 
sented by double decreasing exponential laws so that the 
mass density decreases moving away from the equatorial 
plane and the galactic center or both. 

Taking a system of polar coordinates with origin at 
the galactic center [r, 6, <j>], the height above the equato- 
rial plane is z — r cos 6 and the distance from the galactic 
center along the equatorial plane is R = rsin#, where 8 
is the angle between the polar vector r and the z-axis 
perpendicular to the galactic plane passing through the 
center. The azimuthal symmetry rules out the angle <f) 
and therefore the density laws for the three components 



Pi = poi exp 



r sin 9 



exp 



r Icos 0\ 



where V can be " * " , "M" or "MC" that is stars, dif- 
fuse ISM and star forming MCs. R* d , Rf c , and Rf axe 
the radial scale lengths of stars, MCs and ISM, respec- 
tively, whereas z^, z^ c , z^ 1 are the corresponding scale 
heights above the equatorial plane. To a first approxi- 
mation, we assume the same distributions for stars and 
MCs so that R* d = Rf c and z* d = zf c thus reducing the 
number of scale parameters. Distributions with different 
scale parameters for the three components are, however, 
possible. 

The constants poi vary with the time step. Let us 
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indicate now with tG the age of the galaxy to be mod- 
elled. For the gaseous components only the normalization 
constants pomc (tG) and pom (to) are required as both 
loose memory of their past history and what we need is 
only the amount and chemical composition of gas at the 
present time. This is not the case for the stellar compo- 
nent for which po* (t) is needed at each time < t < tG- 
In other words to calculate the stellar emission we need 
to properly build the mix of stellar populations of any 
age t' = tG — t as result of the history of star formation. 

The normalization constants are derived by integrat- 
ing the density laws over the volume and by imposing the 
integrals to equal the mass obtained from the chemical 
model. In general, the mass of each component M; is 
given by 



Mi 



Poi exp 



r sin 8 

"~rT~ 



exp 



r | cos 8\ 



dV (8) 



The mass of stars born at the time t is given by ^(t), and 
therefore po* (t) will be obtained by using M* (t) = 
Mm (tG) is the result of gas accretion, star formation and 
gas restitution by dying stars. The current total mass 
Mmc (to) is a fraction of Mm (tG)- The double integral 
(in r and 8) is numerically solved for poi to be used in 
eqn.Q 

The galaxy radius R ga i is left as a free parameter of 
the model, thus allowing for systems of different sizes and 
distributions of the components. 

There is a final technical point to examine, i.e. how 
to subdivide the whole volume of a disk galaxy into a 
number of sub-volumes so that the energy source inside 
each of these can be approximated to a point source lo- 
cated in their centers. This requires that the coordinates 
[r, 8, 4>] are divided in suitable intervals. As far as the 
radial coordinate is concerned, test experiments carried 
out in advance have indicated that subdividing the galaxy 
radius in a number of intervals n r going from 40 to 60 
would meet the condition, secure the overall energy bal- 
ance among the sub- volumes, speed up the computational 
time and yield numerically accurate results. 

The number of radial intervals is derived by imposing 
that the mass density among two adjacent sub-volumes 
scales by a fixed ratio pj/pj+i = C, where £ is a constant. 
Upon simple manipulations the above relation becomes 
rj+i — Tj + C- Therefore, the radial grid is equa lly spaced 
in constant steps given by R ga i/n r (Silva 1999). 

The grid for the angular coordinate 8 is chosen in 
such a way that spacing gets thinner approaching the 
equatorial plane, thus allowing for more sub-volumes in 
regions of higher mass density. We split the angle 8 going 
from to n in ng sub- values. We need an odd value for 
77,9 so tha t we have (n g — 1) /2 sub-angles per quadrant. 
Following ISilval lll999f) . the angular distance «i between 
the two adjacent values of the angular grid is chosen in 
such a way that R ga i subtends a fraction / < 1 of the disk 
scale height (za)\ in other words «i = arcsin (fzd/R ga i)- 
Logarithmic angular steps are then imposed A log 8 — 
(2/ (ng — 3)) log (tt /2a\) where ng is determined by the 
condition that the second angular bin near the equatorial 



plan is gai, with g < 3. This implies ng 



21og(ir/2ai) 



+ 3. 



ably finely spaced near = and to get progressively 
broader and broader moving away clockwise and coun- 
terclockwise from 4> — 0. The angular steps are spaced 
on the logarithm scale. As a matter of fact, thanks to the 
azimuthal symmetry it is sufficient to calculate the radia- 
tion field impinging on the volume V (ri, 8i,<j>i = 0) from 
all other volumes V (rt,8k, 4>k)- A grid thinner at <j> ~ 
than elsewhere guarantees the approximation of point- 
like energy sources at the center of the volume elements 
and the conservation of the total energy in turn. 



3.3.2 Early-type galaxies 

The luminosity distribution of early-type galaxies can 
be described by the density profiles o f Hubble, de Vau- 
couleurs and King dKormendvl I1977T) . the most popu- 
lar of which is the Kin g law that yields a finite cen- 
tral density of mass dFroehlichlll982l) .' However, following 
iFioc fc Rocca-Volmerangel 1119971) . we use a density pro- 
file slightly different from the King law in order to secure 
a smooth behavior at the galaxy radius R ga i- We suppose 
that the mass density profiles for stars, MCs, and diffuse 
ISM are represented by the laws 



pi = P0i 



1+ - 



(9) 



where again "i" can be "*", "M" or "MC" (stars, diffuse 
ISM and MCs) and r*, r^ IC , r^ 4 are the core radii of the 
distributions of stars, MCs, and diffuse ISM, respectively. 
The density profile has to be truncated at the galactic 
radius R ga i , which is a free parameter of the model. This 
would prevent the mass M (r) — > oo for r — + oo. In anal- 
ogy to what already made for disk galaxies, the constants 
Po* (t), Pomc (tG) and p 0M (tc) can be found by integrat- 
ing the density law over the volume and by equating this 
value of the mass to the correspondent one derived from 
the global chemical model 



Poi = 



M, 



4tt 



Rgal/rl 



(10) 



(1+X 2 )'' 



dx 



where x = 
as in Sect 



r/rl while poi and Mi have the same meaning 



. The integral is numerically evaluated 
and solved for poi. 

Like in the case of disk galaxies, the last step is to fix 
the spacing of the coordinate grid [r, 8, 4>]. The problem 
will be simpler and the calculations will be faster than in 
the previous case because of the spherical symmetry. The 
spacing of the radial grid is made keeping in mind the en- 
ergy conservation constrain. Also in this case we take a 
sufficiently large number of grid points (n r ~ 40 — 60). 
The condition on the density ratio between adj acent vol- 
umes, pj/pj+i = £ with £ constant iSilvalll999h . implies 



1 + 



R g al 



i I n r 



(11) 



The grid for the angular coordinate (f> is chosen to be suit- 



where usually r c — r*. The coordinate 8 is subdivided 
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into an equally spaced grid, with ng elements in total, 
and 9i — 0, 9 ng — n. Since the distribution is spheri- 
cally symmetric, i.e. independent from 9, we do not need 
a thinner grid toward the equatorial plane. For the az- 
imuthal coordinate <& we adopt the same grid we have 
presented for disk galaxies. 



3.3.3 Intermediate-type galaxies 

Intermediate-type galaxies are characterized by the rela- 
tive proportions of their bulge and disk: they go from the 
early SO and Sa where the bulge is big, to the late Sc and 
Sd where the bulge is small or negligible. In our models, 
we take all this into account by means of different SFHs 
for the disk and the bulge. Furthermore, both in the bulge 
and disk we consider the three components (ISM, MCs 
and stars) in a realistic way. In analogy to what already 
made for purely disk galaxies we adopt a system of po- 
lar coordinates with origin at the galactic center [r, 6,<f>]. 
Azimuthal symmetry rules out the coordinate <j)- 

In the disk, the density profiles for the three compo- 
nents are the double decreasing exponential laws of eqn. 
. Accordingly, we introduce the scale lengths B , 
R¥,b, R¥,b ') z d,B, z^b and Zd*B , where the suffix B indi- 
cates that now the scale lengths are referred to the disk- 
bulge composite model. In the bulge the three compo- 
nents are distributed according to the King-like profiles 
defined in eqn. @, where now the core radii r* B , rff B 
and r^g are referred to the bulge. The constants of nor- 
malization are derived in the same way described in Sects. 
EHm and l3~3~2l The two SFHs of disk and bulge are as- 
sumed to be independent and are simply obtained from 
the chemical models where the mass of disk and bulge are 
specified. The content in stars, MCs and ISM in a given 
point of the galaxy will be simply given by the sum of 
the contributions for the disk and bulge. 

Owing to the composite shape of the galaxy (a sphere 
plus a disk), we have to define a new mixed grid sharing 
the properties of both components, i.e. those described in 
Sect. l3.3~T1 and Sect. |3~1Q1 Let us indicate with Rb the 
bulge radius and with R ga i the galaxy radius. The radial 
grid is defined in the following way. We build two grids 
of radial coordinates, r B ,i and rn,i, one for the disk and 
one for the bulge, in analogy to what we did in Sects. 
13.3.11 and 13.3.21 As the grid of the bulge is not equally 
spaced, but thicker toward the center of the galaxy, we 
take the coordinates r^s of the bulge grid if rt < Rb, 
while if Ti > Rb, we take the coordinates of the disk m,i 
until Rgai- For the angular coordinate 9 we proceed in 
the same way. We build 8b, i and 9o,% as in Sects. 17.3.11 
an d l3~3~31 In this case the disk grid is thinner toward the 
equatorial plane of the galaxy whereas the bulge grid is 
equally spaced, so we take the coordinates 9o,t of the 
disk as long as 9 D ,i+i — 9 D ,i < 9 B ,i+i - B ,i and 9b, i 
elsewhere. For the azimuthal grid there is no problem as 
it is chosen in the same way both for the disk and the 
bulge. 



3.4 The elemental volume grid 

Assigned the geometrical shape of the galaxies, the den- 
sity distributions of the three main components, and 
the coordinate grid (r, 9, <j)) (The number of grid points 
for the three coordinates is n r + l,n$ + l,n^,), the 
galaxy is subdivided into (n r ,ne,n<s>) small volumes V, 
each one identified by the coordinates of the center 
(riv,0jv,<t>kv) given by the mid-point of the coordi- 
nate grid r iV = (n + n+i) /2, 9 jv = {9 J +9 J+1 )/2 
and 4>kv = {4>k + 0fc+i) /2. Thereinafter the volume 
V (riv, 9jv, ®kv) will be simply indicated as V(i,j,k). 
The mass of stars, MCs, and diffuse ISM, in each vol- 
ume are easily derived from the corresponding density 
laws pi (i, j, k) V (i,j, k) where i stands for stars, MCs, 
and ISM. By doing this, we neglect all local gradients in 
ISM and MCs (gradients inside each elemental volume). 
Since the elemental volumes have been chosen sufficiently 
small, the approximation is fairly reasonable. 



4 EXTINCTION AND EMISSION OF THE 
DIFFUSE ISM 

IPiovan et alj fcOOd) presented a detailed study of the ex- 
tinction and emission properties of dusty ISMs. They take 
into account three dust components, i.e. graphite, sili- 
cates and PAHs and the final global agreement reached 
between theory and the ISM extinction and emission data 
of MW, LMC and SMC has been very good. As we are 
now going to include this dusty ISM model in our galax- 
ies, it is wise to briefly summarize here the basic quanti- 
ties and relationships in usage. 

First of all, the total cross section of scattering, ab- 
sorption and extinction is given by 

°v (A) = ira 2 Q p (a, A) — ^^-da (12) 

Ja min , i n H da 

where the index p stands for absorption (abs), scatter- 
ing (sea), total extinction (ext), the index i identifies 
the type of grains, a m i n ,i and a max ,i are the lower and 
upper limits of the size distribution for the i-type of 
grain, uh is the number density of H atoms, Q p (a, A) are 
the dimension-less absorption and scattering coefficients 
dPraine fc Led 11984 lLaor fc Dramel Il993t iLi fc Draind 
l200^1) ~and, finally dnj(a)/da is the dist ribution law of 
the grains (IWeingartner fc DraindfeoOlaf) . With the aid 
of the above cross sections it is possible to calculate the 
optical depth r p (A) along a given path 

t p (A) = g v (A) J n H dl = g v (A) x N H (13) 

where L is the optical path and all other symbols have 
their usual meaning. In this expression for t p (A) we have 
implicitly assumed that the cross sections remain con- 
stant along the optical path. 

Let us name j s x mal1 , j b x i3 and jf AH the contributions 
to the emission by small grains, big grains and PAHs, re- 
spectively. H ow these quan t ities a re calculated is widely 
described in IPiovan et alJ i2006t) to whom the reader 
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should refer for more details. Let us summarize here just 
the key relationships in usage. 

The contribution to the emission by very small grains 
of graphite and silicates is 



■ small 

3x 



aQabs (a,X)B x (T(o)) X 



dP(a) dT 1 rfn(a) 
dT riH da 



da 



(14) 



where dP (a) /dT is the distribution temperature from 
Tmin to Tmax attained by grains with generic dimen- 
sion a under an incident radiation field and B\ (T (a)) 
is the Planck function. Q a bs (a, A) are the absorption co - 
efncients, dn (a) /da is the lWeineartner fc Draina ll200Tal ) 
distribution law for the dimensions, a/; u is the upper 
limit for thermally fluctuating grains, a m i n is the lower 
limit of the distribution. 

The emission by big grains of graphite and silicates 
is evaluated assuming that they behave like black bodies 
in equilibrium with the radiation field. Therefore we have 



• bin 
J\ = T 



a Q ab s (a,X)B x (7») 



1 dn (a) 
riH da 



da 



(15) 



where a max is the upper limit of the distribution and the 
meaning of the other symbols is the same as in eqn. 1141 
The emission by PAHs is given by 



■ PAH 



untie j\ . 

,/PAiT 



/ ( A ) A f"" M0 X 

« da 



a 2 ^Qibs A " [a,X J Sion {^X , A, a) \i+ (16) 
a, A ) Sneu (A , A, a) (1 — xO dadX 



a bs 

*.NPAH 
i abs 



where the ionization of PAHs is taken into ac- 
count l|Weingartner fc Drainel l2001rJl and \i ~ Xi ( a ) 
is the fraction of ionized PAHs with dimension a. 
Sion(x,X,o^ and Sneu (x , A, a^j give the energy 
emitted at wavelength A by a molecule of dimension a, as 
a consequence of absorbing a single photon with energy 
hc/X . a l p"% H and a p *?g are the limits of the distribution 

and I I A ) is the incident radiation field. 



with the radiation field. Two simplifying hypotheses are 
worth being made here: 

(i) The dust contained in a generic volume V does 
not contribute to the radiation field impinging on the vol- 
ume V' . This approximation stands on the notion that, 
in first approximation, owing to the low optical depths of 
the diffuse ISM in the MIR/FIR, dust cannot effectively 
absorb the radiation it emits, except for high density re- 
gions such as MCs, for which dust self-absorption has al- 
ready been taken into account. In other words, the dust 
contained in V' will be transparent to the IR radiation 
emitted by the dust contained in the volume V . There- 
fore, only stars and MCs will contribute to the incoming 
radiation. 

(ii) Following lRvbicki fc Lightmanl dl979t) . the radia- 
tive transfer from a generic volume V to V' is simply cal- 
culated by means of the so-called effective optical depth 
defined by 



T eff = \/ T abs {Tabs + T 3 . 



(17) 



The above relation stands on the following considera- 
tions: scattering increases the absorption, however as 
photons are not destroyed, the effective optical depth 
must be lower than the sum of the optical depths by scat- 
tering an d absorption but higher tha n the one by sole ab- 
sorption dRvbicki fc Lightmanl Il979t) . Although this ex- 
pression for r e ff has been derived for photons randomly 
travelling across an absorbing diffusive medium, so that 
it would strictly apply only to an infinite, homogeneous, 
isotropically scattering medium, we make use of it here. 
The results from the above approximation have been 
compa red bvJSib/a^^^ ( 1998) with those bv lWitt et alj 
l| 1992ft and lFerrara et alJ Jl999ft using Monte-Carlo meth- 
ods to solve radiativ e transfer p r oblem s. T he results fully 
agree with those bv lWitt et alJ il992Tl and lFerrara et all 
( 1999) in the case of spherical symmetry and partially 
disagree in the case of disk galaxies. For these latter the 
quality of agreement dep ends on view ang le between the 
galaxy and the observer dSilva et al.lll99aTl . 
The total radiation field incident on V' is 



n r rifj n% 



H^v) = EEE 

i=l j = l k—l 



V-[j MC (X,V) + 



r 2 (V, V) 



-r (A,n) 



■8) 



= [-Ta//(A,V,T/')] 



5 THE GALAXY SED 

Given the main components of a galaxy, their spatial dis- 
tribution, the coordinate system, and the grid of elemen- 
tal volumes, to proceed further one has to model the in- 
teraction among stars, dusty ISM and MCs to simulate 
the total SED emerging from the galaxy. 

Let us consider a generic volume V' = V(i',j',k') 
of the galaxy: it will receive the radiation coming from 
all other elemental volumes V — V k). The radiation 
travelling from one volume to another interacts with the 
ISM comprised between any two generic volumes. There- 
fore one has to take into account the energy that is both 
absorbed and emitted by the ISM under the interaction 



where the summations are carried over the whole ranges 
of i,j,k with i / / j' and k k'; r 2 (V,V) is the 
value averaged over the volume of the square of the dis- 
tance between the volumes V and V'; r e ff (A, V, V') is 
the effective optical depth from V to V' defined by eqn. 
1171 . which by means of eqn. becomes 



e //(A, V, V') = sj (Jabs (A) X {a abs (A) + a sea (A)) X 

r v(i',j',k') 

n H {l)dl (19) 



V(i,j,k) 



This integral represents the number of H atoms contained 
in the cylinder between V and V . All details on the 



© 2005 RAS, MNRAS 000, 000-000 



Galaxies spectra with interstellar dust 9 



derivation of this quantity and r 2 (V, V') are given in Ap- 
pendix A and B. 

The two terms j MC (A, V) and j* (A, V) are the emis- 
sion by MCs and stars per unit volume of V (i, j, k). They 
are evaluated at the center of the volume element. 

Let us now define two kinds of SSPs: those that are 
already free of the parental cloud and are indicated as 
sspf (classical free SSPs), and those that are still em- 
bedded in their parental dusty molecular clouds and are 
indicated as ssp d (dusty SSPs). 

Let us denote with fd the fraction of the SSP lumi- 
nosity that is reprocessed by dust and with to the time 
scale for this to occur, in such a way that 

C 1 t^t 
fd=< 2 - t/t to < t < 2t (20) 

[ t^t 

where to is a suitable parameter determining the evapora- 
tion time of the parental MCs. Accordingly, the fraction 
of SSP luminosity that escapes without interacting with 
dust is // = 1 - f d . 

The parameter to will likely depend on the properties of 
the ISM and type of galaxy in turn. Plausibly, to will be 
of the order of the lifetime of massive stars. It will be 
discussed in more detail in Sect. 15. ll 

The monochromatic luminosity of a free SSP of given 
age t' — to — t, born at t, and metallicity Z for unit of 
SSP mass is therefore 

rM v 

/ $(M)fx(M, T ',Z)dM 
Li (r, Z) = JM « (21) 

/ $ (M) dM 

Jm l 

Knowing the the monochromatic luminosity of the 
naked SSPs L{ (t', Z), the monochromatic luminosity of 
the dust ens hrouded SSPs Lf (t , Z) has been derived as 
described in lPiovan et alJ (I2006T) . The emission of stars 
and MCs per unit volume, j MC (A, V) and j* (A, V) re- 
spectively, are given by 



we can obtain the emission per unit volume from the 
dusty ISM. At this point the azimuthal and spherical 
symmetries of the galaxy models become very important. 
The emission per unit volume from the dusty ISM cal- 
culated at the center of the volume elements is the same 
everywhere, independently of the coordinate Therefore 
it is sufficient to calculate the dust emission at $ = for 
all possible values of r and 9 on this "galaxy slice". It is 
obvious that the symmetry will be lost when considering 
the total emission from a certain volume element because 
owing to the adopted spacing of the galaxy the elemental 
volumes are not equal. However, as long as we refer to the 
emission per unit volume, the symmetry properties above 
allows us to avoid tedious and lengthy numerical calcu- 
lations. The total radiation field for unit volume emitted 
by a single element is 

jTOT (A) v) = .MC (A) v) +f v) +j ISM (A; v) (24) 

where j ISM (A, V) is the radiation outgoing from a unit 
volume of the dusty diffuse ISM and is given by the sum 
of the contributions from silicates, graphite and PAHs 
described by eqns. (1141 . (1151 . and 11611 . The total out- 
going emission from the volume V is therefore given by 
jTOT ^ x y k v i ous iy different from volume to vol- 
ume. 

The radiation emitted by each elemental volume 
(n r ,ng,n$) has to travel across a certain volume of the 
galaxy itself before reaching the edge, escaping from the 
galaxy, and being detected by an external observer. While 
finding its way out, the radiation is absorbed and diffused 
by the ISM. The external observer will see the galaxy 
along a direction fixed by the angle O, where O = 
means that the galaxy is seen face-on, whereas B = 7r/2 
means that the galaxy is seen edge-on. To this aim, we 
need to determine the properties of the cylinder of mat- 
ter from the center of each element V to the edge of 
the galaxy, along the direction O in order to calculate 
the amount of radiation emitted by the galaxy along this 
direction. Therefore, the monochromatic luminosity mea- 
sured by an external observer is 



j*(X,V)= f G P*(t)L[ (r',Z)dt + 

J2t 

+J t 2t0 (r ~ 1 )p*^ L i( T '' z ) dt ( 22 ) 

and 

j MC (\,V)= f° p.(t)Li(r',Z)dt + 
Jo 

+ J' t0 (2-£)p.(t)L d x (T',Z)dt (23) 

It is worth noticing that luminosity of the MCs de- 
pends only on the luminosity of the young embedded 
stars and not on the mass of molecular gas enclosed in 
the MCs. The factors (2 - t/t ) and (1 - (2 - t/t )) = 
(t/to — 1) follows from relations 12011 and the definition 
of//- 

Once calculated the incident radiation field J (A, V') 



L{X,e) = ^Y.T.zZ v xj TOT (A,V) 

i=l j — l k=l 

xe [-r. // (A,v,e)] (25) 

where T e // (A, V, 0) is the effective optical depth between 
V k) and the galactic edge along the direction 0. The 
details on the derivation of the effective optical depth 
r e ff (A, V, 0) are described in Appendix C. 

5.1 SSPs in usage and parameters of the galaxy 
model 

In this section first we shortly present the libraries of stel- 
lar models, isochrones, and stellar spectra that are used 
to calculate the SEDs of SSPs and galaxies. Second we 
summarize the res ults for SSP i ntrinsically affected by 
dust calculated bv lPiovan et alJ fcOOl l200rj) and in the 
case of young SSPs we describe the effect of the parame- 
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ters of the model. Finally, we present the other parame- 
ters governing the the ISM and chemical model and the 
galaxy geometry. 

(1) Libraries of stellar models and stel- 
lar spectra: W e adopt the set of isochrones 
by iTantalo et alJ il998T) (anticipated in the 
data base for galaxy evolution models by 
iLgj Hierer. Alloin. von Alvensleben. Gallager. Huchra fc et 
(1996)). The underlying stellar models are those of the 
Padova Library calculated wi th convective o vershootin g 
and are amply described by iFaeotto et ail l)l994alltM . 
so that no detail is given here. 

The library of stellar spectra is from iLeieune et alJ 

( 1998), which stands on the Kurucz (1995) release of the- 
oretical spectra, however with seve r al imp ortant imple- 
mentations (see IPiovan et, alJ 120031 . l200fil . for more de- 
tails) . 

(2) Dust enshrouded SSPs As already men- 
tioned, young stars and AGB stars are both surrounded 
by their own dust. The young stars because they are im- 
mersed in MCs, the AGB stars because they eject dust 
on their own. In both cases the UV-optical light emit- 
ted by the stars is absorbed and re-emitted in the FIR. 
But for old ellipticals, in which star formation stopped 
at early epochs because of the galactic wind so that no 
young stellar populations embedded into parental MCs 
are present, for all the other morphological types the im- 
pact of young dusty populations on the galaxy SED has 
to be considered. This indeed is stronger than the one 
caused by dusty AGB stars (see below). The reason of it 
can be reduced to the fact that while for young SSPs all 
the emitted energy interacts with the dusty environment, 
in the case of old SSPs only AGB stars interact with the 
dust. However, as the ranges of wavelength interested in 
the case of young SSPs and dusty AGBs are different, 
the accurate description of the region around 1 — 2fim 
requires that the dusty shells around AGB stars are fully 
taken into account. SSPs for A GB stars and young stars 
with this effect built in are by IPiovan et all (120031) and 
IPiovan et aD (|200fih . 

Intermediate age and old SSPs: In this paper we de- 
cided not to use the SSPs of IPiovan et~aT1 i2003t) where 
the effects of circumstellar dusty shells around the AGB 
stars (of intermediate/old ages) are included. The main 
reason is that dust enshrouded AGB stars are available 
only for three metallicities. Work is in progress to build 
a more complete set of SEDs in which the effects of the 
dusty shells around AGB stars are included. 

Young SSPs: In general dust around young stars 
shifts the UV/optical radiatio n into the M IR/NIR. How- 
ever, the SED of young SSPs IPiovan et al.ll2006^l is gov- 
erned by a number of parameters describing the MC it- 
self, chief among which are the optical depth ry in the 
V band, the scale radius of the cloud R, the C abun- 
dance be in the two log-normal populations of very small 
carbonaceous grains according to the distribution law of 
IWeingartner fc Draind l2001al) , and finally the ionization 
state of PAHs. In addition to these, there is another pa- 
rameter related to young MCs, namely the evaporation 
time to. Let us examine their effect in some detail, 
(i) The MC optical depth ry in the V band. Up to now 
in our library of SSPs only two values are considered, i.e. 



tv = 5 and 35, corresponding to intermediate and high 
optical depths. The effect of changing the optical depth 
is simple: the higher the optical depth, the higher is the 
amount of energy shifted toward longer wavelengths. A 
remark is worth here. Using a certain set of SSPs (i.e. 
a set with given parameters) we implicitly fix the op- 
tical depth. Clearly, the optical depth of MCs not only 
all increases with the cloud mass, decreases with the cloud 
size, and in general increases with the cloud density, but 
also increases with the metallicity. Therefore the ideal 
situation would be the one in which MCs in galaxy mod- 
els cover the whole range of masses, radii, densities and 
metallicities. 

(ii) The scale radius of the MC. Two values are al- 
lowed, R = 1 and 5. They correspond to MCs of different 
compactness. Keeping constant all other parameters, the 
peak of emission shifts t o longer wavelength s at increas- 
ing scale radius (see also lTakagi et al1 l2003). 

(iii) The carbon abundance be per H nucleus fixes the 
abundance of this element in the two log-normal popu- 
lations of very small grai ns in thelWe ingartner fc Dr aind 
l2001al) law. As shown bv lWeingartner fc Draind 1200 1 a ). 
the sole extinction curve is unable to constrain be- It 
provides only an upper limit for be, which at given ra- 
tio Rv = A (V) /E (B — V) is reached when the very 
small carbonaceous grains and PAHs are able to account 
for the ultraviolet bump of the extinction curve. Values 
lower than the upper limit are possible. Without other 
constraints, for instance the IR emission of the region of 
interest, it is not possible to determine be, which plays 
an important role in the IR emission of galaxies. Specif- 
ically, it drives the emission of PAHs in the MIR gener- 
ated by young MCs. The effect of varying be is simple: for 
higher values of be, the PAH emission in the MIR reaches 
higher flux levels. In contrast, the UV-optical flux and the 
shape of the absorbed stellar emission do not depend on 
be, because the extinction curve remains the same for 
different values of be- As the global abundance of C is 
fixed, the distribution of the grains has to compensate for 
the small number of small grains with an higher number 
of big grains. Therefore, for low values of be we expect 
that the emission in the FIR slightly increases (the total 
energy budget has clearly to be conserved). 

(iv) The ionization state of PAHs. The optical proper- 
ties we ha ve adopted are diffe rent for ionized and neutral 
PAHs (see lLi fc Draindl200ll for more details). The ion- 
ization state of PAHs affects the PAH emission profile in 
the MIR. We considered three ionization mod els. The first 
is the one bv lWeingartner fc Drainel l2001bf) . The second 
model adopts for the MC s the same ionization profile 
calculated bv lLi fc Drainel 1200 ll) for the diffuse ISM of 
the MW. The third model simply takes into account only 
neutral PAHs (very low ionization). In other words, we 
have a seq uence of models go ing from ionized to neutral 
PAHs (see IPiovan et al.l l2006. for mo re details). 

(v) to: Following the simple recipe bv lSilva et ail ll998l) . 
this is the scale time required by a new generation of 
stars to g et rid of the parent al MCs (see also lPiovan et all 
l2006h . In lSilva et all ll998l) model, star formation is re- 
duced to a point source scheme, in which stars born at 
the center of dusty MCs slowly e vaporate them to ev en- 
tually shine free. As pointed from lTakagi et ail l2003tl . if 
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the time of the star forming activity is shorter than the 
escaping time scale to , the light coming from young stars 
is complete ly hidden a t the U V-optical wavelenghts. In 
general, in Silva et alJ (|l99Si) model, for times shorter 
than to, no light escapes from dusty star forming re- 
gions and therefore the light observed in the UV-NIR 
should be negligible for any aperture size of the galax- 
ies. This is not what we observe in real galaxies espe- 
cially in their central regions. The discrepancy is cured 
by adopting the more realistic descriptions in which stars 
are randomly distributed inside the MC, and solving the 
radi ation transfer probl em with the ray tracing method 
fsee iPiovan et alJl2006h . We can reproduce in detail in 
this way the inner regions of starbursters like Arp220 
and M82. Even with the new scheme, the parameter to 
drives the amount of energy emitted by the young stars 
that is absorbed and reprocessed by the local dust (i.e. 
in the region immediately around the stars themselves). 
The evaporation is simulated by letting more and more 
energy to escape without being reprocessed by dust. The 
time scale to likely goes from 3 to 100 Myr, i.e. the evo- 
lutionary lifetime of a 100 and 5 Mq respectively. High 
values of to mean that young stars are longer hidden by 
the parental clouds and accordingly much of the light 
they emit in the UV-visible range is shifted to the FIR 
for a long period of time. In such a case a large fraction 
of the IR light emitted by a galaxy could be due to young 
stars still embedded in MCs. The opposite for low values 
of to. 

(3) Parameters of the ISM model are the frac- 
tionary mass of gas in the ISM /m and the metallicity 
Z. 

(a) fivi: The fraction of diffuse gas present in a galaxy 
with respect to the total amount of gas can be inferred 
from observational data on the ratio between the molec- 
ular and atomic gas. In disk galaxies for instance most of 
hydrogen in MCs is in form of H2, whereas the inter-cloud 
medium is mainly made of atomic HI. This parameter 
bears very much on emission and extinction by the dif- 
fuse ISM. The larger /m, the bigger is the amount of 
gas present in the diffuse ISM. Therefore, the emission 
(proportional to the number density of H atoms) will be 
higher. Furthermore, extinction of the diffuse ISM as well 
will increase because of the bigger absorption by dust, 
which is also proportional to the number density of H 
atoms in the diffuse ISM. 

(b) Metallicity Z: The metallicity reached by the ISM 
of the galaxy fixes the properties of the dust, drives the 
intensity of the ISM emission and its ability in extinguish- 
ing the radiation field. We adopted as standard models 
for Z < Zq the three compositions of MW, LMC and 
SMC, that allow us to cover well this range of metallici- 
ties. The dust to gas ratio scales with the metallicity: the 
higher the metal content of a galaxy, the higher is the 
abundance of grains per H atom. This simply implies 
that the ISMs of high metallicity galaxies will both ab- 
sorb and emit more energy. For metallicities higher than 
the solar one (Z > Zq), we keep the relative proportions 
of silicates, graphite and PAHs of the diffuse ISM of the 
MW, scaled to the higher dust to gas ratio. 

(4) Parameters of the chemical models: Chief 
among others are the exponent k and the efficiency v of 



the star formation law, the infall time scale r, and lim- 
ited to the case of elliptical galaxies the parameters fixing 
the gravitational and thermal energy of the gas and driv- 
ing the onset of galactic winds. In the case of galaxies 
composed by disk and bulge, the number of parameters 
increases. There are k B , vb and tb together with ko, 
vo and td, and the bulge to disk mass ratio Mb/Md- 
Finally, we will also consider the case in which galaxies 
may undergo a burst of star formation. The burst is in- 
troduced superimposing it to the star formation history 
law described by eqn. ||SJ. The burst is described by the 
following parameters: the age tsuR at which the burst 
occurs, the enhancement factor vbur amplifying the cur- 
rent efficiency v, and, finally, the duration AtBUR- 

(5) Parameters of the geometrical models: In 
the case of early type galaxies (our spherical case) the pa- 
rameters are r*, , rff , i.e. the core radii for the distri- 
bution of stars, MCs and ISM, respectively. In the case of 
disk galaxies the spatial distributions of the three compo- 
nents are fixed by the radial scale lengths R a , R d IC , and 
R^f , and the scale heights z* d , z d IC , and z d ' , with obvious 
meaning of the various symbols. For composite galaxies 
with disk and bulge the number of parameters increases 
dramatically as we are dealing with two distributions: 
the King law (as for early-type galaxies) for the bulge, 
described with the parameters r * B , r^ B > t^b an d a dou- 
ble decreasing exponential law for the disk, described by 
the parameters R* d B , R% B , Rf )B and z\ B , z d ' B , z d ' B . In 
all types of model, the scale length of the young and old 
stellar populations is the same. Even if this reduces the 
number of parameters, in reality it would be interesting, 
for starburst galaxies in particular, to consider different 
spati al distributions for old and young stars Csee lCalzettil 
2001). Work is in progress to improve upon this type of 
model. The remaining geometrical parameters are: the 
radius of the galaxy R ga i and, if a bulge is present, the 
radius R B of this, and finally the angle of inclination go- 
ing from edge-on to face-on galaxies, i.e. from O = 0° to 
= 90° with respect to the z-axis perpendicular to the 
equatorial plane of the galaxy. 

5.2 Role of the parameters in template model 
galaxies 

In this section we analyze in detail the effect of the various 
parameters on galaxy SEDs. To this aim we calculate 
three template galaxies, i.e. an old early-type galaxy, a 
late-type spiral galaxy, and finally a starburst galaxy. 

5.2.1 Early-type galaxies 

Early-type galaxies are simulated by means of spherical 
models whose total baryonic mass is 10 n Mn and total 
dark matter mass is 5 times higher (see iTantalo et alJ 
1996., for details). For the purposes of this section we 
limit ourselves to consider the typical SFH of an elliptical, 
characterized by a strong initial burst of stellar activity, 
early onset of galactic wind, and quiescence ever since. 
As the model stands on the infall scheme, the SFR starts 
small, increases to a peak value and then declines ever 
since. The time scale at which the peak occurs is of the 
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Figure 1. Basic quantities of the chemical model for a pro- 
totype early-type galaxy as function of the age: the top-left 
panel shows the star formation rate in Mq /yr; the top-right 
panel displays the maximum (Z, solid line) and mean metal- 
licity ((Z), dotted line); the bottom- left panel shows the mass 
of living stars M s t ar (solid line), the gas mass M gas (dotted 
line) , and the total mass of baryons (dashed line) ; finally 
the bottom-right panel displays the ratios Mstar/M^ (solid 
line) and M ga3 /M^ (dotted line). All masses are in units of 
10 11 Mq. Ages are in Gyr. 



order of the infall ti me scale that is fixed at t=0.1 Gyr 
IITantalo et alJll9 96^) which means that the model tends 
to the closed-box approximation but for the occurrence 
of galactic winds. The exponent for the star formation 
law is fe=l. The efficiency of star formation is v = 6. The 
re ason for the ab o ve choice an d more details can be found 
in iTantalo et al.l Jl996lll998l) . The four panels of Fig. □ 
show a few relevant quantities characterizing the galaxy. 

To study the effect of the parameters we consider two 
stages of the history of an early-type galaxy: the age of 
0.15 Gyr, when the SFR reaches the peak and the present 
age of 13 Gyr after billion years of passive evolution. 

In the very early stages of high star formation we 
can assume that dust is mostly concentrated in the dense 
regions of star formation. Following this idea we assign 
a low gas content to the diffuse ISM (/a/ = 0.3, which 
implies that the mass in dusty MCs is about twice mass 
of the diffuse ISM). We adopt t ~ 40 Myr, a long evap- 
oration time, which sounds reasonable for a high density 
star forming environment. Once the galactic wind has 
taken place and star formation has ceased, we can adopt 
/a/ = 1 and to = 0. 

In the 13 Gyr model there is no longer star forma- 
tion and the evolution is merely passive. If the galaxy is 
free of gas and only stars are present, the SED is ex- 
pected to drop off long-ward of about 2\im. However 
elliptical galaxies of the l ocal universe emit in the IR 
iGuhathakurta et alJll986l) . The origin of this flux in the 
MIR/FIR may be due to diffuse dust which emits at those 
wavelengths. Therefore to match this part of the spec- 
trum one has to allow for a small amount of diffuse ISM 
which is likely to exist. By imposing /a/ = 1 for our model 



of age 13 Gyr, we assume that all the gas is in form of 
diffuse ISM. Interesting questions to rise are: how much 
gas can be present today in an elliptical galaxy? What is 
the source of this gas? 

To answer the above questions let us examine how 
the gas content of an elliptical galaxy is expected 
to ev olve with time iGibson fc Matteucc 119971: IChiosi 
2000). According to the classical scenario jLarsonlll974Tr 
all the gas present in elliptical galaxies at the onset of the 
galactic wind is supposed to be expelled from the gravi- 
tational potential well of the galaxy into the intracluster 
medium. Despite the radiative losses, the energy input 
from supernova explosions, stellar winds and dynamical 
interactions overwhelms the gravitational potential and 
makes galactic wind to occur in such a way to happen 
later and later at increasing galaxy mass. Star formation 
suddenly ceases and consequently the energy input by 
Type II (early on) and Type I (afterwards) either stops 
(Type II) or get very small (Type I). All this requires 
about 0.5 Gyr. Subsequently low mass stars (< 2Mq) 
loose mass by stellar wind (during the red giant and the 
asymptotic giant branch phases) thus refueling the galaxy 
of gas in amou nts that are c omparable to those before the 
galactic wind (jchiosi 200oj). What i s the fate of this gas? 

In the most plausible scenario iGibson fc Matteucc! 

1997j) the phase of galactic wind should last for about 
0.5-1 Gyr thanks to the energy input by Type la su- 
pernovae and then stop. The amount of gas lost by low 
mass s tars during t his time interval turns out to be very 
large JChiosilEoobT) and it is not entirely clear how it 
can exceed the gravitational energy. Nevertheless, there 
is general consensus that on a rather short time scale the 
gas may escape the galaxy. Most likely a sort of dynam- 
ical equilibrium is reached in which gas is continuously 
ejected by stars and lost by the galaxy. It may happen 
that a tiny amount of gas is always present in the galaxy, 
thus accounting for the IR emission. 

After the onset of the galactic wind, our chemical 
model is not able to describe this complex situation. Dy- 
ing stars emit lots of gas whose fate is uncertain. Basing 
on the flux level observed in real galaxies, the gas content 
of the galaxy is reduced by a factor of the order of 1000. 

Finally, we are left with the geometrical param- 
eters to fix. For the exponents 7* and jmc a value 
around 1.5 could be taken for elliptic als. The v a lue fo r 
the ISM is more uncert ain. Following iFroe hlichl lll982l) . 
IWitt et alJ Jl992l) . and IWise fc Silval lll996t) we choose 
7a/ — 0.5 — 0.75, that is the ISM is less concentrated 
toward the center than the stars. 

Before galactic wind, at the ages of 0.15, the gas is 
made of molecular clouds and diffuse ISM. Therefore, we 
need the scale lengths of MCs, dusty ISM and stars for 
which we assume r* /C = r* = r^f = 0.5 kpc. Afterwards, 
we need only the scale lengths for ISM and stars. Both are 
kept unchanged (0.5 kpc). In general we will assume r* = 
r^ IC in spheroidals to reduce the number of parameters. 

In Table columns (2) and (3), we summarize the 
set of parameters we have used to model our test ellipti- 
cals. 

At the age of 13 Gyr only the geometrical parameters 
need to be discussed. First we consider the galaxy radius 
Rgai, which is let vary from 20 to 0.5 kpc, i.e. from very 
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Figure 2. Top-left panel: star formation history of the pro- 
totype elliptical galaxy with Ml = 10 11 Mq at the age of 13 
Gyr. Top-right panel: SEDs of elliptical galaxies with radii 
of 20, 10, 5, 3, 1 and 0.5 kpc. The full range of wavelengths is 
represented, from 0.1 to 500^tm. Bottom-left panel: details 
on the UV-optical/NIR region of the SEDs. Bottom-right 
panel: contribution to the total emission of the diffuse ISM 
for the different radii. 



Figure 4. In the top-left and bottom-left panels we show 
the total emission of the prototype elliptical galaxy at the 
age of 13 Gyr. The ratio r*/r^f between the scale lengths is 
kept fixed and two values are considered: r^/r^' 1 =0.1 and 
r */ r c f = 3- In both cases r* is then varied from 0.1 to 9 kpc. 
The corresponding value of r^f is shown. In the right panels 
(top and bottom), the contribution to the total emission of 
the diffuse ISM is shown in detail. 



expanded to very compact systems. In Fig. [5] we show 
the resulting SEDs together with the emission in the UV- 
optical region (bottom-left panel) and the contribution to 
the total flux by the diffuse ISM (bottom-right panel). At 
decreasing size of the galaxy, the optical depth increases: 
the dimension of the galaxy scales linearly, the density of 
matter and the numerical density of H atoms uh increase 
as cc r -3 and, therefore, the optical depth r increases as 
cc r~ 2 . From the top-right and bottom-left panels of Fig. 
121 we can see the effect of it on the UV-optical/NIR part 
of the spectrum: the smaller the dimension, the stronger 
is the extinction of the UV-optical light. For the same 
reason the emission of dust in the FIR becomes stronger 
at decreasing dimensions because of the higher density 
(bottom-right panel of Fig.^l. The emission, in fact, lin- 
early depends on the numerical density of H atoms. 

As the mass of the physical components of the galaxy 
and all their scale lengths are fixed (in our example all 
equal to 0.5 Kpc), in models of smaller size more gas, 
stars and dust are stored in the outer regions with re- 
spect to the larger size models. In such a case the mat- 
ter distribution tends to become uniform, a limit situa- 
tion reached when r c 3> R ga i- For compact objects, the 
peak of emission in the FIR also shifts toward shorter 
wavelengths. This is due to the higher temperature of 
the grains, thanks to the combined effect of the higher 
local emission and the presence of more sources in the 
outer regions (bottom-right panel of Fig.|5|. 

Other geometrical parameters are the scale length 
r* c and the ratio r*/r^' that at fixed r* c gives r^ 1 . We 
calculate a grid of models characterized by seven values 
of r*, i.e. 0.1, 0.5, 1, 3, 5, 7 and 9 kpc {R ga i is fixed at 
10 kpc), and ten values of the ratio r*/r^ , i.e. 0.005, 
0.05, 0.1, 0.2, 0.4, 1, 3, 10, 50 and 100, for a total of 170 



models covering large intervals of r* and r c . The results 
obtained for various combinations of these parameters 
are presented in figs. |3]and2] In Fig. |3] we show the re- 
sults obtained keeping fixed r* a and leaving free the ratio 
r c/ r c 1 > or equivalently . Four values of r* axe consid- 
ered. For each of these on the left panel we show the total 
galaxy flux at varying r^r, whereas on the right panel we 
display the corresponding contribution of the diffuse ISM 
to the total flux. For the two smallest values of r*, 0.5 
and 1 kpc, at varying r^f the emission in the FIR grows, 
reaches a peak value, and then declines. This can be ex- 
plained in the following way. For r ^ <C r * and r£ 1 3> r* , 
the diffuse ISM is concentrated either in the inner or in 
the outer regions of the galaxy. In both cases the spatial 
distributions of the ISM does not favors the interaction 
with the stellar radiation: the FIR flux increases as r;r 
grows from very small values to r^f ~ r*, which repre- 
sents the best condition in which the intensity of the local 
radiation field is strong in high density regions. Moving to 
r^ 1 2> r* the flux decreases, because the dust is concen- 
trated in regions of weaker radiation field. The emission 
peak shifts to longer wavelenghts because for small r*, 
as r^f increases, the dust is confined in regions of lower 
radiation field intensity. For the two highest values of r * , 
5 and 7 kpc, the flux simply increases as dust goes from 
being concentrated in the inner regions of the galaxy, to 
being shifted toward the outer regions with higher den- 
sity of stars. 

In Fig. 0| we keep fixed the ratio r*/r^ and let r* c 
vary from 0.1 to 9 kpc. The result is the same for both 
cases no matter whether r^/r^ 1 < 1 or > 1. The FIR 
emission is stronger when both r* and are small, 
i.e. when stars and dust occupy the same inner region 
of the galaxy. With high values of the two radial scale 
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Table 1. Columns (2) through (5): parameters for a sample of model galaxies, namely two ellipticals (E) at the ages of 13 and 
0.15 Gyr, a disk (D) at the age of 13 Gyr and a starburst (SB) at the age of 13 Gyr. Columns (6) through (11) best-fit parameters 
of the models we have adopted to match the SEDs of the spiral galaxies M100 and NGC 6946, the elliptical galaxies NGC2768 
and NGC4494, and finally the local starburst galaxies Arp220 and M82. 
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4 Age in Gyr. 2 Distance of the galaxy in Mpc. 3 Baryonic mass of the galaxy in units of IO^Mq. 4 Evaporation time in Myr of 
young SSPs from dusty parental regions. In the 13 Gyr models and the galaxies NGC2768 and NGC4494 with no ongoing star 
formation, the parameter is set equal to 0. fractionary gas content in the diffuse ISM. "Ratio between the core radii of stars 
and ISM in the spherical models. 7 Scale radii in kpc of stars and MCs in the spherical models. 8 Radial scale in kpc of the diffuse 
ISM in the disk models. 8 Radial scales in kpc of stars and MCs in the disk models. 9 Vertical scale in kpc of the diffuse ISM in 
the disk models. 10 Vertical scales in kpc of stars and MCs. n Radius of the galaxy in kpc. 12 Infall time scale in Gyr. 13 Efficiency 
of star formation. 14 Age of the beginning of the burst for starburst models in Gyr. 15 Length of the burst for starburst models. 
1B Multiplying factor for the star formation for starburst models. 17 Inclination angle. 
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Figure 3. Left Panels: the SEDs of the test elliptical galaxy at the age of 13 Gyr for r* = 0.5, 1, 5 and 7 kpc. For each of these 
the ratio r*/r^f is varied from 0.005 to 50. The corresponding values of r^f are also shown. Right Panels: the contribution to 
the total emission by the diffuse ISM. The arrows show the evolution of the cirrus emission going from the smaller to the higher 
values of r^f . All the other parameters are kept fixed. 



lengths, stars and dust are more dispersed across the 
galaxy, therefore the FIR emission is weaker and the peak 
shifts toward longer wavelengths. 

To show the effects of to, the time scale for a young 
SSP to evaporate the MC in which it is embedded, and 
f M , the fractionary gas mass in the diffuse ISM, we take 
the 0.15 Gyr old model galaxy, at the peak of star for- 



mation. In Fig. |5] we show the SEDs for several values 
of to. The effect of this parameter is strong and straight- 
forward: the longer this time scale, the higher is the ob- 
scuration of the UV-optical light because the energy of 
newly born young and luminous objects is long trapped 
by the dusty environment (top-right panel). Simultane- 
ously, more energy is emitted in the FIR by the dusty 
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Figure 5. Top-left panel: star formation history of the pro- 
totype elliptical galaxy of Ml = 1O 11 M0 with a mark on the 
modelled age at 0.15 Gyr. Top-right panel: SEDs of ellip- 
ticals at varying the time scale to from 1.5 Myr to 100 Myr. 
The considered values are 1.5, 4.5, 15, 30, 60 and 100 Myr. All 
the other parameters are kept fixed. We also show the SED of 
the same galaxy calculated with the classical EPS technique. 
Bottom-left panel: contribution to the total emission by 
dusty MCs at varying to- Bottom-right panel: contribution 
to the total emission of the diffuse ISM for the same set of 
time scales to- 

MCs regions (bottom-left panel). In any case it is clear 
from the top-right and the bottom-right panel that an 
increase of to, corresponds to a decrease of the emission 
by the diffuse ISM, whose peak also shifts toward longer 
wavelengths. The explanation is that increasing to the 
UV-optical intensity of the local radiation field dramati- 
cally drops, leaving a cooler and weaker mean radiation 
field to heat the dust. 

Finally, we compare in Fig. |H] the SEDs of our test 
galaxy at varying the amount of gas in the diffuse ISM, 
i.e. /m. The remaining gas is stored in MCs. As shown by 
the top panels, increasing the amount of gas in the diffuse 
ISM has little influence on the global SED. The reason is 
that the galaxy is at the peak of star formation, where the 
attenuation and IR emission are clearly dominated by the 
star forming MCs. Similar albeit smaller effects of /m are 
expected in starburst galaxies with high star formation 
rate. In any case, as expected, increasing the amount of 
gas in the diffuse ISM yields stronger attenuation of the 
light emitted by old and young stars (top and bottom left 
panels) . The emission by the diffuse ISM (which anyhow 
scarcely contributes to the total flux) increases with the 
amount of gas and also becomes cooler because of the 
weaker local radiation field heating the grains (bottom 
right panel). 

5.2.2 Late-type galaxies 

We consider a disk galaxy of 10 1 Mq at the age of 13 
Gyr. In our model disk galaxy (i.e. with infall scheme and 
the Schmidt law of star formation) and likely in real disk 



Figure 6. Top panels: SEDs of test ellipticals at varying the 
fraction of gas in the diffuse ISM /m with respect to the total 
amount of gas. The range of values goes from 0.1 to 0.9. All 
the other parameters are kept fixed. The considered values are 

0. 1, 0.3, 0.5, 0.7 and 0.9. We also show the SED of the same 
galaxy calculated with the classical EPS technique. Bottom- 
left panel: contribution to the total emission by dusty MCs at 
varying /m- Bottom-right panel: contribution to the total 
emission of the diffuse ISM for the same set of time scales Jm- 

galaxies as well, the star forming activity never ceases. 
Due to the lack of significant galactic winds in disk galax- 
ies, the gas remain for ever in the disk and continuously 
refuels star formation. This means that our model has to 
take all possible components into account even at the age 
of 13 Gyr: young stars just formed and still embedded in 
their parental molecular cloud, bare stars of any age and 
diffuse ISM. The task is more complicate than with the 
spheroidal models due to the lower degree of symmetry 
(only the azimuthal one). 

The choice of the three main parameters driving the 
star formation and chemical enrichment of a disk galaxy, 

1. e. the infall time scale r, the exponent k, and the ef- 
ficiency v of the star formation rate, rests on the fol- 
lowing considerations . In the star fo rmati on rate, k typ - 
ically varies from 1 (Schmidt 1963|)to 2 (|LarsorJll99lT) . 
In our model for a lat e -type galaxy the choice is guided 
by the lPortinari et al] (Il998f) results for the Solar Vicin- 
ity. Therefore we adopt k = 1, and v = 0.7, whereas 
for the infall time scale we assume r = 4. The resulting 
star formation rate starts small, grows to a maximum at 
about 3.5 Gyr, and then gently declines to the present 
day value. 

The temporal evolution of a few physical quantities 
characterizing the model galaxy is shown in Fig. 

To proceed further other parameters must be spec- 
ified. First, we fix the geometrical parameters, i.e. the 
radial and vertical scale heights, for which we make the 
most simple choice, i.e. R* d = R\ IC = R^f = 5 kpc and 
z d ~ Z ¥ C ~ z ¥ = 0.4 kpc (the spatial scales for the 
three components are the same) in agreement with typ- 
ical value s found for local spir al galaxies such as M100 
and M51 feeckman et alll99d) . Second, we need to spec- 
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Figure 7. Basic quantities of the chemical model for the test 
late-type galaxy as function of the age: the top left panel shows 
the star formation rate in Mq /yr; the top right panel displays 
the maximum (Z, solid line) and mean metallicity ((Z), dot- 
ted line); the bottom left panel shows the mass of living stars 
Mstar (solid line), the gas mass M gaa (dotted line), and the to- 
tal mass of baryons (dashed line) ; finally the bottom right 
panel displays the ratios Mstar/M^ (solid line) and M 3as /M^ 
(dotted line). All masses are in units of 10 11 Mq. Ages are in 
Gyr. 
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Figure 8. Top panels: SEDs of prototype disk galaxies at 
varying the radius of the galaxy R ga i, going from 0.5 to 20 Gyr. 
All the other parameters arc kept fixed. We also show the SED 
of the test galaxy calculated with the classical EPS technique. 
Bottom-left panel: contribution to the total emission com- 
ing from dusty MCs at varying R ga i- Bottom-right panel: 
contribution to the total emission of the diffuse ISM for the 
same set of radii R ga i- 



ify the fraction /m of the gas present in form of diffuse 
ISM. The task is not easy and to a certain extent be- 
yond the aims of this study because it would require a 
multi-component description taking into account all heat- 
ing and cooling processes that transform atomic likely hot 
gas into cool molecular gas. So, as already presented in 
Sect. 15. ll we consider /m as a parameter, and, for the sake 
of illustration, we adopt /m = 0.5. This means that half 
of the gas is in the diffuse medium and half in the MCs. 
Third, for the timescale to we assume ~ 5 Myr which 
roughly corresponds to the lifetime of the most massive 
stars. This means that the heavily obscured part of young 
stars lifetime is rather short. For the sake of illustration 
we consider the case with inclination angle = 60°. 

In Table [T] column (4), we summarize the set of pa- 
rameters we have used to model our test disk galaxy. 

In the same way we did for the ellipticals we calcu- 
lated a sequence of models with R ga i going from 20 to 
0.5 kpc, from very expanded to very compact systems. 
In Fig. |H]we show the resulting SEDs (top panels a and 
0) also with the spectrum calculated with the classical 
EPS. In the bottom panels we show the emission in the 
UV-optical region (bottom-left panel) coming from MCs 
and the contribution to the total flux coming from the 
diffuse ISM (bottom-right panel). The effect of the ra- 
dius is the same as for ellipticals: for smaller galaxies, 
higher optical depths, the stronger will be the extinction 
of the UV-optical light. For the same reason the emission 
of dust in the FIR will be stronger for the smaller dimen- 
sions because of the higher density (bottom-right panel 
of Fig.®. 

In the case of disk galaxies the degree of symmetry is 
lower than for spheroidal systems as only azimuthal sym- 
metry is conserved. In Fig.|§] we show the SED at varying 
the radial scale lengths R$ and R d . In the four left pan- 
els we keep fixed Rff and vary R d , the opposite in the 
four right panels. Examining the four panels on the left, 
we notice that once fixed R d , the attenuation becomes 
weaker going from low to high values of R d (panels a and 
0). This is simply due to the fact that growing the value 
of R a , more stars are distributed in the outer regions of 
the galaxy and thus they are less obscured. It is worth 
noticing how in our model disk galaxy the attenuation is 
partially due to young MCs and partially to the diffuse 
ISM. The ultimate reason of this is the low evaporation 
time for dusty regions we have chosen for the model. In 
panel 7 we see the effect of attenuation on the emission 
of MCs: again it is stronger for R d < R d J , because for 
these values young dusty SSPs are confined in the inner 
region (R d IC = R d ) ■ The effect of varying the star scale 
on the diffuse ISM emission (panel 5) is that this emis- 
sion becomes weaker and cooler at the increasing of R d , 
because more stars are distributed in the outer regions 
of the galaxy weakly heating the ISM. Passing now to 
the four right panels the above effects are reversed. R d is 
fixed and the obscuration of the stellar light by the dif- 
fuse ISM (panels a and 7) grows at the increasing of R d , 
because the ISM shifts toward the outer regions more and 
more wrapping the stellar component. In the same way 
the emission of the ISM becomes stronger at the increas- 
ing of R^ 1 , because dust is more evenly distributed in the 
regions occupied by stars. 
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Figure 9. Four Left Panels: in panels a and /3 we show the SEDs of the 13 Gyr test disk galaxy with fixed R¥ = 5 kpc. The 
radial scale length R* d is varied from 1 to to 9 kpc. The following values are considered, 1, 2, 3, 4, 5, 7 and 9 kpc. The SED obtained 
with the classical EPS technique is also shown in panel a for the sake of comparison. The dashed line in panel /3 represents the 
contribution of MCs to the total flux in the FIR. In panel 7 we show the effect of varying R* d on the MCs flux. In particular the 
uv-optical/NIR spectral region is shown. Finally, in panel 6 we plot the contribution of the diffuse ISM to the total emission. Four 
Right Panels: the same as above but for fixed R* d and varying Rj~- R?f goes from 1 to 9 kpc. 
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Figure 10. Four Left Panel: in panels a and f} we show the SEDs of the test disk galaxy of 13 Gyr for the fixed value of the 
vertical scale length z¥ = 0.4 kpc. The scale length z* d is varied from 0.1 to to 1.5 kpc. The considered values for z* d are 0.1, 
0.25, 0.4, 0.55, 0.7, 1 and 1.5 kpc. The SED obtained with the classical EPS technique is also shown in panel a for the sake of 
comparison. The dashed line in panel [3 represents the emission of the MCs in the FIR,. In panel 7 we show the effect of varying 
z d on the MCs flux. In particular the uv-optical/NIR, spectral region is shown. Finally, in panel 8 we plot the contribution of the 
ISM to the total emission. Four Right Panels: the same as above but for z* d fixed and z^ varying from 0.1 to 1.5 kpc. 



In Fig. 1101 we show the effect of the vertical scale 
lengths. As we did for the radial scales, we consider two 
cases: first we fix z d ' and let z* d vary, second we fix z d 
and let z\ l change. In the four left panels z^ is fixed. At 
increasing z d the extinction becomes lower, because more 
stars are distributed in the outer regions of the galaxy (cv 
panel). This effect of growing attenuation is the same 
on young dusty SED (panel 7) that in this model are 
distributed with the same vertical scale of bare stars. The 



ISM emission tends to be stronger for lower values of z d 
(panel S). In the four right panels of Fig. llUl the vertical 
scale of stars z d is kept fixed. In this case the effect of 
varying the other parameter is smaller: in particular in 
panel 8 we can see as the emission of the diffuse ISM 
grows to a maximum and then decreases. The maximum 
is reached when the scale lengths of stars and dust are 
similar. If z& <C z d , or z d * 3> z d the emission tends to 
be lower. 



© 2005 RAS, MNRAS 000, 000-000 



18 L. Piovan, R. Tantalo & C. Chiosi 




X [fi m] X [[i m] 

Figure 11. Top-left Panel: star formation history of the 
prototype disk galaxy with = 10 11 Mq. The age of the 13 
Gyr model whose SEDs are examined in detail is marked. Top- 
right Panel: SEDs of the galaxy at varying the evaporation 
time scale to from 1.5 Myr to 100 Myr. The considered values 
are 1.5, 4.5, 15, 30, 60 and 100 Myr. All the other parameters 
are kept fixed. We also show the SED calculated with the 
classical EPS technique. It steeply declines for wavelengths 
longer than 2fim. Bottom-left Panel: contribution to the 
total emission from dusty MCs at varying to- Bottom-right 
panel: contribution to the total emission of the diffuse ISM 
at varying tg. 

In Fig. 1111 we show the effect of varying the evap- 
oration time to of young dusty SSPs. The longer the 
time scale, the higher is the effective extinction of the 
UV-optical light because the energy of newly stars is ob- 
scured for a longer time (top-left panel), and the higher 
is the amount of energy shifted toward the FIR by MCs 
(bottom-left panel). As for the elliptical galaxy of 0.15 
Gyr, at increasing to the emission of the diffuse ISM de- 
creases and becomes cooler, because the intensity of the 
local radiation field is weaker. 

Finally, we compare the SEDs of our disk galaxy at 
varying the amount of gas in the diffuse ISM with respect 
to that in dusty MCs, using the parameter /m- This is 
shown in Fig. 1121 The effect of /m can be easily seen in 
the top panels (a and 0). Increasing /m, i.e. the amount 
of gas in the diffuse ISM, strongly affects the SED. The 
higher the amount of gas, the stronger is the obscura- 
tion, and the higher is the flux in the FIR. As expected, 
increasing the amount of gas in the diffuse ISM makes 
stronger also the effect of attenuation on the light from 
young stars (bottom left panel). The emission of the dif- 
fuse ISM becomes also cooler with more gas in the diffuse 
ISM, because of the weaker radiation field heating the 
grains (bottom right panel). 

5.2.3 Starburst galaxies 

Starburst galaxies are objects that show a recent and 
transient increase in SFR by a large factor (ten at least). 
The burst is often confined to a few hundred parsecs near 
the nucleus, although bursts extending to wider regions 
are easy to find. The high SFR of starburst galaxies is 
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Figure 12. Top panels: SEDs of model disk galaxies at vary- 
ing the fraction of gas in the diffuse ISM f^i respect to the 
total amount of gas. The range of values goes from 0.1 to 0.9. 
All the other parameters are kept fixed. The considered val- 
ues are 0.1, 0.3, 0.5, 0.7 and 0.9. We also show the SED of the 
prototype spiral calculated with the classical EPS technique. 
Bottom-left panel: contribution to the total emission com- 
ing from dusty MCs at varying Jm- Bottom-right panel: 
contribution to the total emission of the diffuse ISM for the 
same set of /m- 



of great interest, because it is the local analog of high 
redshift galaxies during their formation, involving strong 
star formation in dust rich environments. 

The light emitted by newly born stars in these star 
forming galaxies is absorbed by dust and re-emitted in 
the FIR/MIR. One of the main problem is to disentan- 
gle whether the dominant energy source heating dust is 
a starburst or an AGN. In our models, AGNs as energy 
source is not yet included. Future work is planned to in- 
clude the AGN contribution thus providing a more real- 
istic description of these systems. In Tabled column (5), 
we summarize the set of parameters we have adopted to 
model our prototype starburst galaxy. 

To simplify things we take a spherical model for 
which we adopt a long infall time scale and relatively 
low star formation efficiency, i.e. t — 9 and v = 1. In 
this model star formation is maximum at about 3 Gyr, 
slowly declines up the the present time (13 Gyr) and 
never ceases. It is a sort of late type galaxy but for the 
spherical symmetry. As already mentioned the specific 
shape is irrelevant for the purposes of this experiment. 
At the age of 12.95 Gyr (not long ago) we introduce a 
short burst whose intensity is 30 times stronger than the 
current SFR at the age of 12.95 Gyr. 

Star formation in starbursters is tightly related to 
reddening and obscuration which render the detailed in- 
terpretation of their continuum and emission-lines very 
complicate. Starbursters are known to possess a flat 
obscuration law (w ith no bump) as pointed out by 
ICalzetti et al.l lll994T) . who considered a mix of stars and 
dust and took the effects by scattering into account. This 
obscuration law is very puzzling, because starbursters are 
typically objects of high metallicity and so we would ex- 
pect an extinction law similar to that of the MW. It 
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is worth noticing, however, that the dust obscuration 
in galaxies is not strictly equivalent to dust extinction 
in stars. The latter measures the optical depth of the 
dust between the observer and the star, whereas the for- 
mer corresponds to a more general attenuation, to which 
many effects may concur such as those by extinction, 
scattering, and geo metrical distrib ution of the dust rela- 
tive to the emitters l|Calzetti200ll f . Because of this, there 
is nowadays much debate o n the origin of the fla t law. Is it 
due to geometrical effects l|Granato et al.ll200(t ). or pecu- 
liar distribution and the composition of the dust in star 
forming environments? Furthermore, there is some evi- 
dence th at reddening by gas and stars are systematically 
different (Calzctti 200L), as stars are on average less red- 
dened than ionized gas. Another point is that both stars 
and ionized gas are contrasted by a foreground-like dis- 
tribution of dust. This could be due to dust being pref- 
erentially associated to the star forming regions. In ad- 
dition to this, the evidence of different amounts of gas 
mass along different lines of sight, could mirror a den- 
sity and redde ning structure. These considerations led 
ICalzettil l|200iri to propose a model aimed at explaining 
many observational constraints: newly born stars form in 
a central high density region immersed in a bath of older 
stars and dust. 

In our models we can easily deal with this bath of 
old stars and dust: setting r* = (in spherical models) 
and/or R* d = R d ' (in disks), the region of dust coincides 
with the region in which stars are present. We can also 
allow for different scales between old and newly born stars 
so that centrally confined burst can be superimposed to 
a "normal" galaxy, setting r* / rf c or R* d / R d ' c ' . 
However, as in all this paper, we will simply put r* = 
r* /c or R* d = R d 4C , leaving this point aside because it 
would require a deep investigation beyond the purposes 
of this study. 

As for our starburster galaxy we have adopted the 
spherical symmetry, all scale lengths are the same as in 
model for early- type galaxies. In Fig. 1131 we show the ef- 
fect on the total SED of varying the galaxy radius R ga i- 
The effect is the same as already found for the model 
elliptical and disk galaxies: the smaller the galaxy, the 
higher is the optical depth and the stronger is the atten- 
uation of the UV-optical light. The emission of dust in 
the FIR is stronger for compact objects because of the 
higher density (bottom-right panel of Fig. 1131 

In Fig. 1141 we show how the total SED changes at 
varying the scale radii. Their effects on the emission of 
diffuse ISM and young dusty MCs are also highlighted. 
Two cases are considered. First, the ratio r'/r^ 1 is fixed 
and r* is let vary over a wide range of values (four left 
panels). Second, r* is fixed and is let change (four 
right panels). In both groups of four, panels a and f3 dis- 
play the total SED, the old SED without dust, the contri- 
bution of MCs for the extreme values of the parameters, 
whereas panels 7 and 5 show the contribution of MCs 
and diffuse ISM to the total flux. In the four left panels 
of Fig. lMl we notice that the attenuation is stronger when 
r* a and rjr are the lowest, both for stars (panel a) and 
MCs (panel 7). The emission of diffuse ISM grows for 
more compact systems, because the radiation field heat- 
ing the grains is stronger in the central regions of high 
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Figure 13. Top panels: SEDs of the model star-burst galaxy 
at varying R ga i from 0.5 to 20 Gyr. The values are the same 
as for the model elliptical and spiral. All the other parame- 
ters are kept fixed. The SED derived with the classical EPS 
is also shown. Bottom-left panel: contribution to the to- 
tal emission by dusty MCs at varying R ga l- Bottom-right 
panel: contribution to the total emission by the diffuse ISM 
at varying R gal . 
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Figure 15. Top panels: SEDs of the model star-burst galaxy 
at varying the fraction of gas in the diffuse ISM /m with re- 
spect to the total amount of gas. All the other parameters 
are kept fixed. The considered values for /m are 0.1, 0.3, 0.5, 
0.7 and 0.9. We also show the SED calculated with the classi- 
cal EPS technique. Bottom-left panel: contribution to the 
total emission by dusty MCs at varying Jm- Bottom-right 
panel: contribution to the total emission by the diffuse ISM 
at varying f M . 

density (panel 8). In the four right panels where r* is 
fixed, there is a small effect of varying r c , because the 
emission is clearly dominated from the MCs component. 
For higher values of r^ 1 , the system is more obscured 
(panel 7) , because there is more dust in the outer regions 
to screen the stellar light. 

We compare the SEDs of our model starburster at 
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Figure 14. Four Left Panels: in this group panels a and fi show the SEDs of the model starburst at the age of 13 Gyr for fixed 

value of the ratio r* /r^' 1 =0.2. The radial scale length r* is let varied from 0.1 to to 9 kpc. The extreme values are marked with 

the corresponding j™ ■ The SED obtained with the classical EPS is also shown for the sake of comparison. The contribution of 

MCs to the total flux is also shown for the extreme values. In panel 7 we show the MC flux in the UV-optical/NIR and in panel 

8 we plot the contribution of the diffuse ISM to the total emission. Four Right Panels: the same as above but for r* fixed and 
M 

varying r * . 
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Figure 16. Top-left panel: star formation rate of the model 
star-burst galaxy of Ml = lO 11 A/0 on which the age of 13 Gyr 
is marked. Top-right panel: SEDs at varying the evaporation 
time scale to, namely 1.5, 4.5, 15, 30, 60 and 100 Myr. All the 
other parameters are kept fixed. We also show the classical 
SED (dotted line). Bottom-left panel: contribution to the 
total emission by dusty MCs at varying to- Bottom-right 
panel: contribution to the total emission by the diffuse ISM 
at varying to. 

varying the amount of gas in the diffuse ISM. The re- 
maining part of the gas is distributed among the dusty 
MCs. In Fig. ll5l we compare SEDs obtained with different 
amounts of gas in the diffuse ISM, varying the parame- 
ter /m- Increasing the amount of gas in the diffuse ISM 
makes stronger the extinction of the light from old and 
young stars, whereas it makes higher and cooler the FIR 
flux, because of the higher density and the weaker radia- 
tion field. 



Finally, in Fig. EH we show the effect of varying the 
evaporation time scale of young dusty SSPs. The effect is 
the same as in our model elliptical of 0.15 Gyr (see Fig.|SJ 
For longer time scales, the attenuation of the UV-optical 
light becomes stronger (top-right panel), and the amount 
of energy shifted toward the MIR/FIR by MCs becomes 
higher (bottom- left panel). 



6 LATE-TYPE GALAXIES OF THE LOCAL 
UNIVERSE 

In this section, using our models we seek to reproduce 
the SEDs of two late-type galaxies of the Local Universe, 
namely M100 and NGC6946. 

In Table Q columns (6) and (7), we summarize the 
parameters characterizing the models that best match 
the properties of the two galaxies under consideration. 
Part of the parameters are based on observational hints, 
the others are suitably varied to get agreement between 
observational and theoretical SEDs. Specifically, the ge- 
ometrical parameters and distances are from literature 
and kept constant. The gas mass in the ISM and MCs 
is let vary within the range indicated by current obser- 
vational data. The star formation efficiency v and infall 
time scale r are also let vary around the typical values 
currently estimated for spiral galaxies. Finally the evap- 
oration timescale to and the set of dusty SSPs in our 
library are free parameters. 

M100. This Sbc spiral galaxy is one of the most 
important members of the Virgo Cluster, characterized 
by two huge and luminous spiral arms and many other 
smaller ones. The redshift of the galaxy, taken from NED 1 

http://nedwww.ipac. caltech. edu/. 
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Figure 17. SED of the modelled Sbc spiral galaxy M100 
at the age of 13 Gyr (continuous line). We represented 
also the old SED with classical EPS (dotted line), the 
emission of the diffuse ISM (dot-dashed line) and the 
emission from young dusty SSPs (dash ed lin e). Data for 
M100 are taken from iBuat et alj Il989l UV). iDonas et alj 
119871 UV), RC3 catalogue of | d e~yaucouleurs et^ldT " i^}92l 
UBVRI), Ide Jong fc van der Kruitl il994l. BVRIHK) 
2MASS J.Iarrett et alJ 120031 JHK), I RAS JSoifer et alj 
119891: iMoshir fc et alj ll99Cft and, finally, IStark et al.l Jl989t 



FIR). 



(Nasa/Ipac Extragalactic Database) is z = 0.00524. Us- 
ing the Hubble constant Ho = 7 2 km/s/Mpc, i t corre - 
sponds to a distance of 21.8 Mdc. Ishaolev et, alJ l|200ll ). 
who performed a thorough search in literature of the 
distances for the spiral galaxies of their sample, re- 
port for NGC 4321 the distance of 16.1 M pc, based on 
Cepheid distance bv lFreedman et alJ Jl994l) and H = 75 
km/s/Mpc. We adopt here for the distance the average 
of the two estimates, i.e. 19 Mpc. The major and minor 
diameters of the galaxy in arcmin (NED) are 7.4' x 6.3'. 
Taking the average dimension of 7' , we obtain th e linea r 
radius R aa i = 19 K p c. Fo llowing |£hieiMn_^l_alJ (|l 



iKnapen fc Beckmanl dl996h and iBeckman etafl dl996 

we adopt the inclination angle O = 27° and we assume 
that both stars and dust have the same radial and vertical 
scale lengths of 5 and 0.5 kpc, respectively. 

The star formation history of M100 is derived as- 
suming the total baryonic mass of 2 • 1Q 1x Mq, v = 0.7 
and r = 4 Gyr. With this choice, the star formation 
never stops: it starts small, grows to a maximum and 
then smoothly declines. We remain with the following 
parameters to adjust: the gas mass in the diffuse ISM 
and in MCs, the evaporation time to and the library of 
young dusty SSPs. We may also slightly vary the radial 
and vertical scale lengths with respect to current esti- 
mates. For t h e gas masses we adopt the values given by 



Young^al 



Young et al 



1989) based on HI and H2, According to 
19891) . Mr 2 = 0.77 • M g , where M g is the 
total gas mass. It means that this galaxy is dominated 
by the molecular component. A good fit of the obser- 
vational SED is obtained adopting the fractionary gas 



mass in MCs /m = 0.35. For the evaporation time we 
find to — 5Myr, consistent with a dust-poor star form- 
ing medium in which the first supernovae explosions in 
the newly born stellar populations evaporate the parental 
gas. The library of young SSPs best suited to M100 (and 
the other spiral as well) is characterized by the optical 
depth t — 35, R — 5, and high b c abundance. Nothing 
can be said for the ionization state of PAHs. We adopt 
the case with the complete treatment of the ionization 
state. Only the detailed spectrum in the MIR region may 
help clarifying the issue. It is worth noticing that for spi- 
rals, more extended and cooler MCs with R = 5 fit the 
FIR emission better than MCs with R = 1. In Fig. El we 
show our best fit of the observational data. The result is 
remarkably good. 

NGC 6946 is a Sc/Scd nearby galaxy, highly ob- 
scured by the interstellar matter of our galaxy, as it is 
quite close to the Galactic plane and it is seen nearly face- 
on. However, the i nclination angle is uncertain due to the 
global asymmetry llBlais-Quelle^te et al.l2Q Q4|). We adopt 
= 35° (see also iBo nnarel et al.lll988l : ICarignan et ail 
Il990l : l0ev fc Kennicuttll990l) . The distance to NGC 6 946 
is uncer tain and going from 5 (|de Vaucouletirslll97flh to 
10 Mpc jRogstad et alll973l). Our d i stance to this galaxy 
has been taken from TshapIev et alJ l|200lf) who give 5.5 
Mpc as a mean value based on previous studies. From 
the NED online catalogue, the diameters of the galaxy 
in arcmin are 11.5' x 9.8'. Adopting the average dimen- 
sion of 10 .5', we obtain the radius R ga i = 10 Kpc. Ac- 
cording to lTacconi fc Yound jfiT)Mi i. the galaxy has some 
evidence of a spiral structure extending beyond 20 kpc 
with HI emission out to 30 kpc. However they adopted 
the distance of 10.1 Mpc, th at is two times l onger than 
the recent value proposed bv IShaolev et all (120011 ) that 
we have adopted here. Taking 5. 5 Mpc, the radius is 
Rgai = 13 Kpc consistent with the lShaplev et all <l200ll) 
results and the NED diameters. The r adial scale lengths 
of sta rs and gas are also taken from iTacconi fc Younel 
( 1986), however rescaled to our shorter value for the dis- 
tance to the galaxy. We adopt the radial scale length of 
5 kpc, and for the sc ale height we use the same value 
of 1 kpc proposed bv ISilva et alJ il998l) . As far as the 
fraction of total gas embedded in young MCs i s con- 
cerned, the question is controversial. lYoung et all (Il989t> 
reports Mhi/Mh 2 = 2.18: this means that only 1/3 of 
the gas is in molecular form, in goo d agreement with 
the results bv lYoung fc Knezekl il989t) on the molecular 
to atomic gas rati o for the morphological types Sc/Scd. 
iTuffs et all Jl996l) found that the bulk of the FIR lu- 
mino sity arises from a diffuse disk component. How- 
ever, iDevereux fc Younsl lll993l) reports NGC 6946 as a 
galaxy where the molecular gas dominates the interstel- 
lar medium and the thermal emission is explained mainly 
with a warm dust compone nt heated from young m assive 
stars, a resu lt confirmed i n [Mal hotra et aD (|1996). A re- 
cent work bv lWalsh et alJ <l2002l) confirms the importance 
of the molecular component in NGC 6946, almost as mas- 
sive as the atomic one, giving the ratio Mh 2 /Mhi = 0.57, 
which is high with respect to other galaxies of the same 
morphological type. We find that the observational SED 
is best reproduced if the molecular component is about 
as massive as the atomic one. The star formation history 
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Figure 18. SED of the modelled Scd spiral galaxy NGC6946 
at the age of 13 Gyr (continuous line). We represented also 
the old SED with classical EPS (dotted line), the emission 
of the diffuse ISM (dot-dashed line) and the emission from 
young dusty SSPs ( dashed line). Data for NGC 6946 are taken 
from iRifatto et alJ il99l UV 1. Ide Vaucouleurs et all jl99l 
BV), ISO iRoussel et aljboOll MIR ). IR AS iRice et alJ ll98Sl 
MIR and FI R1. lEngargiolal il99ll FIR'). Ibevereux fc Yound 
Il993l. FIR^.lTuffs et alJ il996ft : iTuffs fc Gabriel <2003l FIIVl 
and lSilval il999ft . 

of NGC 6946, that exhibits a moderate starburst activ- 
ity, has been calculated with a baryonic mass similar to 
M100 (1.2 • 10 n M Q , v = 0.7 and r = 5). 

In Fig. 1181 we show the result of our fit to the data 
for NGC 6946. The evaporation time to is very short, 3 
Myr, and the library of dusty SSPs is the same as for 
M100. As compared to M100, the agreement with the 
MIR is better, even if MIR spectra would allow a better 
comparison between theory and observations. 



7 EARLY- TYPE GALAXIES OF THE 
LOCAL UNIVERSE 

In this section we present two old early-type galaxies of 
the Local Universe, namely NGC 2768 and NGC 4491, 
and our attempts to reproduce their SEDs. In Table Q 
columns (8) and (9), we summarize the values of the pa- 
rameters we have chosen. The number of parameters is 
smaller than for disk galaxies. First because of the higher 
degree of symmetry, and second because parameters like 
to and /m have no role. The star formation ceased long 
ago after the onset of the galactic wind, no young stars 
in dusty MCs are present (to = 0) and all gas (if any) is 
in the diffuse ISM (/a/ = 1). There is however, a new pa- 
rameter to consider, i.e. the fraction f r of the gas contin- 
uously ejected by dying and evolved stars (mostly RHB) 
that is still retained in the galaxy (see Sect. l5.2~T1 for de- 
tails). 

NGC 2768 is an elliptical gala xy of mo rphological 
type E6 at th e distance of 21 5 Mp c llTakagi et alJl2003l) . 
As noticed bv lTakaei et alJ i2003l) . this galaxy has been 
detected in the sub-mm range of wavelengths thus pro- 
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Figure 19. SED of the modelled E6 elliptical galaxy 
NGC2768 at the age of 13 Gyr (continuous line). We rep- 
resented also the old SED with classical EPS (dot-dashed 
line) and the emission of the diffuse ISM (dashed line) . 
Data for NGC 2768 are taken fr om lLongo et all Il99ll 
far-UV). jde Vaucouleurs et al.lll992l UBV). NED dat abase, 
iFrogel et alJ Il97gl. JHK1. IRAS iMoshir fc et alJ Il990l. MIR 
and F IR) and NED database and, finally. IWiklind^^Ienkel 
Il995l. sub-mm). 



viding a better constraint on the dust emission. Data for 
this galaxy have been taken from the literature as listed 
in Fig. EH As in iTakaei et alJ (I2003I) . all the data have 
been corrected where it was necessary to be consistent 
with the galaxy observed as a whole, even if there is a 
certai n amount of error for the IUE data of lLongo et alJ 
il99lT) that correspond to a small aperture in the cen- 
ter of the galaxy. The diameters of the galaxy in arcmin 
(NED) are 8.1' x 4.3' and using an average dimension of 
6', we obtain a radius of the galaxy of about R ga i = 20 
Kpc. 

The parameters chosen to model the SFH are as fol- 
lows: the infall time scale is r=0.1 Gyr, the exponent for 
the s tar formation law is fc=l as in lTantalo et alj il996L 
Il99ct) . whereas the efficiency is v—2. In Fig. 1201 we show 
the SFH with other quantities derived by the chemical 
code. The onset of the galactic wind takes place at about 
1.1 Gyr. The scale lengths of stars and gas have been 
fixed as for the prototype elliptical of Sect. I5.2~T1 

Finally, in Fig. EU we show the theoretical SED and 
its comparison with the observational data. Agreement is 
excellent from the UV to FIR. Note the strong reduction 
factor for the gas content. 

NGC 4494 is a roughly spherical elliptical galaxy 
of the morphological type El. For this galaxy the NED 
catalog reports the redshift z — 0.00451 corresponding to 
a distance of abo ut 19 kpc, using Ho = 72 km/s/Mpc. 
iTemi et all i2004l) report from the LEDA catalog a dis- 
tance of about 21.28 Mpc. We choose the average distance 
between the two determinations, i.e. 20 Mpc. The diam- 
eters of the galaxy in arcmin (NED) are 4.8' x 3.5' and 
using an average dimension of 4.2', we obtain a dimension 
of about R g ai = 12 Kpc. 
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Figure 20. Basic quantities of the chemical model for the el- 
liptical galaxy NGC2768 as function of the age: the top left 
panel shows the star formation rate in Mq/jt; the top right 
panel displays the maximum (Z, solid line) and mean metal- 
licity {(Z))\ the bottom left panel shows the mass of living 
stars Mstar (solid line), the gas mass M gas (dotted line), and 
the total mass of baryons (dashed line); finally the bot- 
tom right panel displays the ratios M s t ar /Mr, (solid line) and 
Mqcs/Ml (dotted line). All masses are in units of 10 11 A/q. 
Ages are in Gyr. 
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Figure 21. SED of the modelled E6 elliptical galaxy 
NGC4494 at the age of 13 Gyr (continuous line). We rep- 
resented also the old SED with classical EPS (dot-dashed 
line) and the emission of the diffuse ISM (dashed line) . 
Data for NGC4494 have been taken f rom lRifatto et alj Il995t 
far-UV). Ide Vaucouleurs et alj Il992l. UBV), N ED database 
Gavazzi fc Bosellil 1 19961. UBVJHK1. 2M ASS Ijarrett et alj 
20031 JHK). IRAS jMoshir fc et al"lll99l MIR and FIR) and 
ISO iTemi et alj|2004h . 



The SFH and the evolution in metallicity of 
NGC4494 are calculated in the same way as for 
NGC2768, only with a slightly different value of the bary- 
onic mass. The scale lengths are the same as for the other 
elliptical. In Fig. 1211 we show the result of our fit. The 
agreement is good. The main reason of uncertainty is the 



correc tion we made fo r the IUE UV - data of lRifatto et alj 
<1995ft . As pointed in iTakaei etafl (l2003ft . only knowing 
the true UV profile it is possible to properly correct these 
IUE data that cover only a small region in the central re- 
gion of the galaxy. Finally, the same remark on the gas 
content made for the other elliptical can be made also 
here. 



8 STARBURST GALAXIES 

Finally, we present models for two well known and thor- 
oughly studied star-burst galaxies of the local universe, 
namely Arp220 and M82. In Table [TJ columns (10) and 
(11), we summarize the values of the parameters we chose 
to reproduce the SEDs of these galaxies. We adopt the 
spherical geometry to describe both objects. As for ellipti- 
cals, the number of parameters to deal with is much lower 
than for disk galaxies. The situation, however, is different 
from the case of ellipticals because now the parameters 
/m and particular to play in the fit of the observational 
data. 

Arp220. Arp 220 is the brightest object in the 
Local Universe. There is nowadays the general con- 
sensus that Arp 220 is a starburst-dominated galaxy 



Genzcl et al. 


1998; 


Lutz et al. 


1999 


iRiffODOulou et alj 


ll999l:lTran et alJl200ll). Recentlv. ISooon et alj (120041) re- 



analyzing the ISO MIR spectrum of Arp 220 confirm the 
starburst-dominated hypothesis, suggesting that the IR 
luminosity should be be probably powered by the star- 
burst activity in extremely dense regions, even if the AGN 
contribution cannot be definitely ruled out because of the 
high extinction. 

The redshift of the galaxy, taken from the NED 
database is z — 0.01813. Using the Hubble constant 
Ho — 72 km/s/Mpc, it corresponds to a distance of 
about 76 Mpc. Th is value is fully con sistent with the dis- 
tances proposed in lSoifer et alj lll987l) and in lSpoon et alj 
<l2004t) . The diameters of the galaxy in arcmin (NED) 
are 1.5' x 1.2'. Adopting the mean value of 1.35', the 
radius of the galaxy is of about R ga i = 16 — 17 Kpc. 
IWvnn- Williams fc Becklinl il993tl showed that almost all 
the MIR flux of Arp220 comes from a small central re- 
gion of about 5 aperture and this c oncentration of the 
MIR emission has been confirmed bv lSoifer et al] <ll999l) 
comparing the fluxes at a fixed MIR wavelength and vary- 
ing the beam. For this reason we can quite safely use the 
MIR data from small apertures even modelling the whole 
galaxy. 

The mass of molecular gas H2 in Arp220 according 



m ArpZZU according 
Solomon et alJl997l: 



to CO estimates is a bout 3 TO Mp, 

IScoville et al"lll997l) . lMundell et ail feOOll) estimated the 
nuclear column densities of H2 to be of the order of 
Nh 2 = 2 — 4 ■ 10 22 cm -2 . These values are comparable 
to the mean HI column densities that are of the or- 
der of Nh = 1.5 ■ 10 20 T S , with the spin temperature be- 
tween 100 and 200 K. These column densities are likely 
low er limits because of the uncertainty in the values of 
T a IlKulkarni fc Heiles 1988), in the abundance of CO 
(IFrerking et al] |l982). in the excitation conditions and 
optical depth, which all suggest us to adopt similar con- 
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Figure 22. SED of the model for the starburst galaxy 
Arp220 at the age of 13 Gyr (continuous line). We repre- 
sent also the old SED with classical EPS (dot-dashed line), 
the emission of the diffuse ISM (dotted line) and, finally, 
the emission of dusty MCs (dashed li ne). Optical and near 
IR d ata are taken from RC3 cat a log |d c Vaucou leurs et alJ 
Il99l UBV) and Ijarrett et all fcOOl 2MASS - JHK). 
The MIR fl uxes are fro m Jspinoglio et ail <199El MIR), 
ISmith et alJ <1989t MIR). iKlaas et ail Il997t fr om MIR to 
FIR) IWvnn- Williams fc Becklinl Jl993t MIR). iTran et alJ 
l200ll. 5 -16 ^m ISOCAM-CVF sp ectrum). FIR and radio data 
are from lRigopoulou et_ai[ jl996j. s ub-mm ). lEales et ail I1989L 
UKIRT sub-mm], iDunne et all fcOOOt SCUBA sub-mm), 
iDunne fc EaleJ fcOOll. SCUBA sub-m m). Other MIR/FIR 
data are taken from Tspoon et ail 120041) . 



tents for both atomic and molecular gas. We fix /m = 0.5. 
The SFH history of Arp220 is modelled adding to the 
current SFH (which peaked in past and ever since de- 
creased) a very strong and short burst. The burst is ob- 
tained by increasing the star formation of a factor of 60. 
In Fig. I22l we show our best fit of Arp220, the agreement 
is good. However, there are two points that need to be 
clarified. First, we have no data of Arp220 in the near and 
far UV. As outlined b v lTakaei et alJ J2003l) . the data by 
iGoldader et alJ (120021) are puzzling because, owing to the 
large angular size of Arp 220 filling the field of view of the 
instrument, the flux level of the sky (to be subtracted) is 
highly uncertain. Second, to get satisfactory agreement 
with observational data we had to extend the parame- 
ter space, namely the scale radius has been lowered to 



R = 0.5 and the t v =40, 



ISturm et ail (Il996l) and|Ge 



this one in agreement with 
;el et ail ( 



jenzel et ail J1998I) . The ioniza- 
tion of PAHs is fully considered, but it is worth noticing 
that the ionization model does not bear very much on 
the extinction curve we have adopted thanks to the low 
carbon abundance and low contribution of PAHs. 

There is another interesting point to note. Arp220 is 
modelled as a 13 Gyr old galaxy with a strong burst su- 
perimposed to a spiral-like SFH that reached a maximum 
in the past and gently declines. This long history of star 
formation slowly increases the metallicity well above the 
solar value. For these high metallicities values, one should 
(has to) use a MW-like extinction curve for dense regions 



Figure 23. SED of the model for the starburst galaxy M82 
at the age of 13 Gyr (continuous line). We represent also 
the old SED with classical EPS (dot-dashed line), the emis- 
sion of the diffuse ISM (dotted line) an d the emission of the 
dusty MCs (d ashed line). Dat a are from|d^V i aucouleurset^iJ 
119921. UBV). I Johnson! <196fl VRIJ KL) Jjarrett et alJ <2003l 
2MAS S - JH K). ISoifer et all Il987t I RAS'). iGolombek et ail 
119881. IRASI. lRice et al.Hl98sl IRASI.lKlein et al.Hl988l FIR 
and sub-mm). and lF^rst^r^chrelbTret^lTl2^03Ti 



to calculate the SED of young SSPs emerging from dusty 
regions. This type of extinction curve is characterized by 
a high abundance of carbon. We tried to fit the obser- 
vational SED of Arp220 using the lowest available value 
of b c , but the MIR emission was always too high and 
the 10/im absorption features of silicates not as deep as 
required. To obtain a good fit we had to use for young 
dusty SSPs the extinction curve of the S MC, which is 
poor of carbon. Similar result was found bv lTakaei et alJ 
(2003). There seems to be a point of contradiction, be- 
cause the carbon-poor extinction curve of the SMC is 
characteristic of a low metallicity environment, whereas 
here we are dealing with a high metallicity one. Although 
this point deserves careful future investigation, a plausi- 
ble explanation could be that the high rate of type II su- 
pernovae due to the strong burst alters the composition of 
dust, growing the amount of silicates with respect to that 
of carbonaceous grains. Another hypothesis is that this 
high energy output of the star forming process destroys 
molecules and small grains thus altering the distribution 
of the grains. Finally the geometry of the system could 
play a role. 

M82. The central region of the irregular starburst 
galaxy M82 seems to have suffer ed a strong gravitational 
inter action about 10 s years ago jForster Schreiber et alJ 
120031) with its companion M81 and it shows a re- 
markable burst of star formation activity. M82 is still 
strongly interacting with Ad 81 as pro ved by the com- 
mon envelope of n eutral hydrogen jYun et ail 119931 : 
llchikawa et alJ Il995t). The classical distanc e to M82 is 
3.25 ±0.20 Mpc jTammann fc Sandagdll96g|) that makes 
of M82 the nearest and most studied starburst galaxy 
(see IShopbell fc Bland-Hawthornl Il998i and references 
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therein). For M82 data for the emission of PAHs in 
the MIR region a re available. The ISO-SWS data by 
ISturm et alJ 112000. private communication) cover a wide 
region going from 2.4 to 45/im with different parts of a 
SWS full grating scan that are observed with different 
aperture sizes, going from 14 x 20 to 20 x 33 . The 
problem is that with these small apertures it is not pos- 
sible to observe the whole galaxy, but only part of it. In 
our case the south-weste r n star formation lobe. Recently, 
IForster Schreiber et alJ j2QQ3» P r i va t e communication) 
observed M82 in the MIR range between 5.0 and 16pm 
with ISOCAM on board of ISO. The total field of view 
is 96 x 96 . It covers almost entirely the MIR sources 
of M82 and it is more suitable to be used to model M82 
as a whole. The mass of molecular gas for M82 is of the 
order of 1 — 2 ■ 10 8 Mp, , as derived from CO determinations 
(lLo et al.lll987t IWild et alJll992l : IForster Schreiber et alJ 
2001). The esti mates of the total ma ss of gas are of the 
order of 10 9 M Q jSolinger et al.lll97/H . For this reason we 
use a parameter /m ~ 0.8 — 0.9, with almost all the gas in 
the diffuse ISM. In Fig. I23l we show the best fit obtained 
for M82. The agreement with the data is good and the 
main features of PAHs are well reproduced. The young 
dusty SSPs used for this fit are characterized by the fol- 
lowing parameters: R = 1, tv = 8, b c = 2 ■ 10~ 5 and MW 
average ionization of PAHs. An interesting point is that 
we found it difficult to reach the low level of flux in the 
UV region. Even if data are highly uncertain, this could 
be due to the mod el of MCs with unif orm distribution of 
dust, 7] = 1000 fsee lPiovan et alJ^O Oel. we have adopted. 
It might be possible that for M82, MCs with a shell-like 
distribution of dust with more dust in the outer regions 
and therefore with stronger attenuation of the flux from 
young stars, are better suited to obtain a good fit in that 
spectral region. 



9 DISCUSSION AND CONCLUSIONS 

In this paper, improving upon the standard EPS tech- 
nique, we have developed theoretical SEDs of galaxies, 
whose morphology goes from disk to spherical structures, 
in presence of dust in the ISM. Properly accounting for 
the effects of dust on the SED of a galaxy increases the 
complexity of the problem with respect to the standard 
EPS theory because it is necessary to consider the distri- 
bution of the energy sources (the stars) inside the ISM 
absorbing and re-emitting the stellar flux. This means 
that the geometry and morphological type of the galaxy 
become important and unavoidable ingredients of the 
whole problem, together with the transfer of radiation 
from one region to another. The emergent SEDs of our 
model galaxies have been shown to reproduce very well, 
even in details, the observational data for a few test galax- 
ies of different morphological type. The model is versatile 
and applicable to a large range of objects of astrophysical 
interest at varying the star formation and chemical en- 
richment histories, the geometrical shape or morphology 
of the galaxies and the amounts of gas and dust contained 
in their ISM. 

Before concluding, it is worth mentioning a few 
points of weakness that could be improved . First, the 
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chemi c al mo dels we have adopted are from lTantalo et alJ 
ll996j, Il998|), where as the chemical yields are from 
iPortinari et alJ lll998l) . These models are state-of-the-art 
in the study of the chemical evolution of galaxies. How- 
ever, they do not include a proper description of the for- 
mation / destruction of dust as for instance recently devel- 
oped bv lDwekl |l998, 2005]). Even if the dust content can 
be related to the metallicity of the galaxy (see Sect. 15. ll . 
the relative proportions of the various components of the 
dust would require the detailed study of the evolution of 
the dusty envir onment and t he complete information o n 
the dust yields <Dweldll998L l2005t iGalliano et al]l2005T) . 
This would lead to a better and more physically sounded 
correlation between the composition of dust and the star 
formation and chemical enrichment history of the galaxy 
itself. All this is missing in most galaxy models in which 
dust is considered. The problem may be particularly se- 
vere for high metallicity environments. 

Second, the models for disk galaxies with central 
bulge need to be tested against SEDs of local galax- 
ies of intermediate type going from SO to Sb/Sc, trying 
to match some observa tional constra i nts, l ike the UV- 
optical average colours feuzzon 3 12001120051) . 

Finally, many other physical ingredients can be im- 
proved and/or considered. Just to mention a few, the 
inclusion of the recent models of thermally pulsing AGB 
sta rs with varyin g molecular opacities in the outer lay- 
ers [Marig^ 1120021) . the extension of the SED to the radio 
range, and the simulation of the nebular emission. 
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APPENDIX A: MATTER IN THE CYLINDER 
BETWEEN V AND V 

To calculate the number of H atoms contained in the 
cylinder of matter between two generic volumes V = 
V {i,j, k) and V = V k') one has to know the total 

gas mass in the cylinder 

f V> 

/ PM (1) dl (Al) 
Jv 

where pu is the density profile of the diffuse medium, 
given by the eqns. J7| or @ or both, according to whether 
we are dealing with a disk galaxy, a spherical galaxy or 
a disk plus bulge one. 

Normalizing the density to the central value, i.e. 
p' = pm/pom, and the radial distance to the radial 
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scale length for disks or core radius for spheroidal galax- 
ies as appropriate, i.e. x = r/R^ 1 or x = r/ty, 
we may express the density profiles by means of di- 
mensionless quantities. For a disk galaxy we have 
p'{x) — exp (—xsin 0) exp (—a; |cos#| ■ ii^/z^) which 
depends only on the ratio R$ and zji , whereas for a 
spheroidal galaxy we have p'(x) = [l + (a;) 2 ] " /M ■ 
In this way it is possible to calculate once for all the 
mass contained in the cylinders between any two generic 
volume elements of the galaxy, independently from the 
mass of the galaxy, given the coordinate grid (nr, n9, n$) 
and ratio R^f / (for disk galaxies only) or exponent 7m 
(for spheroidal galaxies). The integral of eqn. IAl| is then 
numerically evaluated. 



APPENDIX B: THE DISTANCE R 2 

Following [Silval (Il999l) the distance r 2 {i,j,k,i' ,f ,k') 
shortly indicated as r 2 (V, V), which is de- 
fined as volume averaged value of the square 
of the distance between V(i',j',k') and all the 
points belonging to V(i,j,k), is r 2 (V, V) = 

d 2 (V, V) r 2 sin 9d9d$>dr/V where d 2 (V, V) 

is given by d 2 (V, V ) = (ry sin 6^ — r sin 9 cos $) 2 + 
(r sin 9 sin $) 2 + (ry cos 6^/ — rcosO) 2 and the volume V 



by V = J J J v r 2 sin 9d9d$dr. 



APPENDIX C: MATTER IN THE CYLINDER 
BETWEEN V AND THE GALACTIC EDGE 

The matter contained in the cylinder between a generic 
element V and the edge of the galaxy extinguishes the 
radiation emerging from the volume element, travel- 
ling across the galaxy up to its edge, from which it 
escapes toward an external observer. To determine this 
mass we proceed as follows. Firstly, since the emission 
jTOT ^ Q £ i [| nas Deen defined per unit 

volume and the problem has azimuthal symmetry, we 
resize the grid of azimuthal coordinates by taking many 
equally spaced bins. Using this resized grid, we introduce 
the Cartesian coordinates (Xv,Yv,Zv) of the center 
of the volume element by means of the transformations 
Xy — Xiv sin#jv cos$ fe v, Yy = xty sin9jy sm$ k y 
and Zy = Xiv cos9jv , where (x^, 9jv, &kv) are the 
polar coordinates of centre of the volume V(i,j,k). 
The starting reference system is in spherical coordinates 
(O, x, 9, $) and corresponds to a Cartesian system 
(0,X, Y, Z). Let us now consider another reference 
system (O' , x' , 9', $') to which another Cartesian system 
(O' , X' ,Y' , Z') would correspond. The transformations 
from one Cartesian system to the other is given by 
X = X' + X ',Y = Y' + Y , and Z = Z' + Z Q ,, 
where (X G i, Y g >, Z G ') are the coordinates of the centre 
O' of the new system in the old one. We take the new 
reference system (O' , X' ,Y' , Z') with the origin O' 
in (Xo>, Y i , Z >) = (0, Xiv sin9jv sin&kv, 0), where 
OO' = Xiv sin 9jv sin & k y. This represents a translation 
of the origin O along the Y-axis to a new origin O' so that 
the centre of the volume lies in the plane (0',X G ', Zq 1 )- 



The Cartesian coordinates of the volume centre in the 
new system will be X' v = Xv, Yy — Yv — Yq' = and 
Z'y = Zv ■ Applying the inverse relationships to pass from 
Cartesian to spherical coordinates, the x' and 9' of the 

volume centre are x' iV = \J (X v ) 2 + (Y v ) 2 + (Z v ) 2 = 
xw \J {sm9j V cos$fcv) 2 + (cos^v-) 2 and 9j V = 
arccos (z v /^J (X v ) 2 + (Y v ) 2 + (Z v ) 2 ^j 



arccos cos 



UjV 



smfjv cos$ fc v) 2 + (cos9jv) 2 



Obviously the new coordinate (p' k is equal to or ir, be- 
cause the volume belong to the plane (O' , Xq 1 , Zq')- For 
x't < (j>' k = 7r, whereas for x\ > cj>' k = 0. To avoid use- 
less complications, the grid of azimuthal coordinates has 
been redefined in such a way that x\ is always different 
from zero. We proceed now to determine the radius of the 
circular section of the galaxy coincident with the plane 
(O' , Xqi , Zqi ) . This will depend on the type of galaxy un- 
der consideration. We get x' G = \J {R ga i/ Rq) 2 — Yq, = 



(Rgal/Ra) - (x iV smOjvsm$kv) ■ 
where R ga i is the galactic radius and Rq is equal to Rj 1 
for disks and to r* f for spherical galaxies. Let us now 
consider an observer looking at the galaxy from the view 
angle Q with respect to the equatorial plane of the galaxy. 
For the sake of simplicity we place the observer on the 
plane (O' , Xq' , Zq')- Therefore, = corresponds to a 
galaxy seen edge-on, whereas O = n/2 to the case face- 
on. Thanks to the azimuthal and equatorial symmetries 
we make take $ — and Z' ^ 0. 

Associated to the volume center with coordinates 
x'j , 9j, <E>^J there will be a point P located on the galaxy 

edge with coordinates P (x'c, 9' G &g) where 9' G and $g 
are still unknown. Let us first calculate &' G . If O — 
all the points at the galactic edge have $q — 0. If 
= 7r/2 two cases are possible: for x\ < =>■ & G = ir, 
whereas for x\ > =>■ &' G = 0. In the general case 
with 0<6<7r/2we have to calculate the equation of 
the straight line with angular coefficient m = tg (it — 0) 

passing through the point Q(X',Z') — Q ^0, — R G ^j ■ 

We obtain Z' = — X'tan© — x' G . We calculate now 
{Z' v + X' v tan 9 + x'q) where {X' v , Z v ) are the coordi- 
nates of the volume center. If {Z' v + X' v tan Q + x' G ) > 
we have & G = 0, whereas if (Zy + X' v tan0 + x' G ) < 
we get <E>' G = 7T. 

The derivation of 9' G is slightly more complicate be- 
cause five cases are possible. In any case, it is a mat- 
ter of lengthy trigonometrical manipulations. The line 
for the center of the plane (O'X'Z') with inclination Q 
has equation Z' — — X'tan0. This leads us to define the 
parameter A to check whether the volume V in the plane 
(O'X'Z') falls above or below the line Z' = -X' tan 6. If 
6 = tt/2 =>• A = 0, whereas if < 6 < n/2 =>■ A = Z'y + 
X'y tanG. Finally, let us introduce the ratio \ — x' iv /x' G . 
The following cases are then possible. If X' v > and A > 
=> then 9' G = tt/2 + 9 - arcsin [xsin (§ - 9 + 9' jv )^ . 

If X'y > and A < => then 9' G = tt/2 + 9 + 
arcsin [xsin (9' jV - tt/2 — 9)] . If X' v > and A = 
then 0'g = tt/2 + G. If X'y < and A > => 
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then e' G = tt/2 + 9 - arcsin [xsin - 6 - Ojv)] ■ 
If X' v < and A < we have three solutions: 6' G = 
tt/2 + 9 + arcsin [xsin (e - tt/2 + 0^)] , 6^ = §tt - 

6 — arcsin j^xsin ^6 — 7r/2 + 9 jy j| and #q = 7r, depend- 
ing on whether + tan 6 + x' G is greater, smaller 
or equal to 0, respectively. 

Once determined the spherical coordinates of 
P {x'giQ'g&g) m the translated system, by means of 
an inverse transformation of coordinates we can ob- 
tain the Cartesian coordinates P (Xp,Yp, Zp) in the 
old system of coordinates: Xp — x' G sin 6' G cos $q , 
Yp = x' G sin 6' G sin $' G + :riv sin cos $ fc y and Zp — 
x' G cos9' G . Having eventually derived the Cartesian coor- 
dinates of the point P (Xp, Yp, Zp) on the galactic edge 
and those of volume center (Xv, Yy, Zy), the calculation 
of the mass in the cylinder comprised between V and 
the galaxy edge is trivial and can be straightforwardly 
performed as described in Appendix A. 
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